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ABSTRACT 


A  complementary  code  pair  has  the  property  that  its  autocorrelation 
function  contains  no  range  sidelobes.   This  property  makes  complementary 
codes  attractive  for  use  in  signal  waveform  design  for  detection  and 
tracking  systems.  An  extensive  list  of  the  properties  of  these  codes 
has  been  formulated  by  M  J.  Golay  and  additions  were  made  by  S.  Jauregui. 
However,  the  ambiguity  functions  of  complementary  codes  have  not  been 
previously  formulated. 

This  thesis  considers  t»he  ambiguity  functions  of  complementary 
codes  for  two  cases: 

1.  Narrow-band  analysis 

2,  Wideband  analysis. 

In  the  narrow-band  analysis  it  is  assumed  that  the  doppler  effect 
does  not  significantly  alter  the  envelope  of  the  transmitted  signal. 
Siebert's  definition  of  the  ambiguity  function  is  utilized  in  the  formu- 
lation of  the  narrow-band  complementary  code  ambiguity  function. 

In  the  wideband  analysis  it  is  shown  that  the  doppler  effect  sig- 
nificantly alters  the  envelope  of  the  transmitted  signal.   This  alteration 
of  the  envelope  of  the  transmitted  signal  becomes  significant  when  the 
product  of  the  two  factors,  the  relative  velocity  between  the  transmitter 
and  the  target,  and  the  time-bandwidth  product,  approaches  one-half  the 
velocity  of  propagation  in  the  medium.   Siebert's  ambiguity  function  is 
extended  to  include  the  wideband  case. 
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CHAPTER  I 
INTRODUCTION 

The  detection  problem  in  modern  search  radar  systems  is  generally 
concerned  with  the  optimization  of  four  major  factors.   These  four 
factors  ares 

1)  Maximization  of  the  effective  range  of  the  system. 

2)  Minimization  of  the  number  of  false  alarms  due  to 
impulse  noise, 

3)  Maximization  of  range  and  velocity  resolution. 

4)  Accurate  measurement  of  the  target's  position  and  velocity 
without  ambiguity. 

The  solution  of  this  optimization  proglem  is  predicated  primarily 
on  a  suitable  choice  of  a  transmitted  signal  waveform.   To  emphasize 
this  point  let  us  analyze  the  four  factors  that  are  involved  in  the 
optimization  problem. 

The  first  factor  can  be  analyzed  from  the  standpoint  of  the  radar 
range  equation.   The  prominent  result  of  this  analysis  is  that  the 
system  designer  may  increase  the  range  of  the  system  by  increasing  the 
duration  of  the  transmitted  pulse. 

The  second  factor  can  be  disposed  of  by  setting  a  high  enough 
threshold  such  that  very  few  noise  impulses  exceed  the  threshold,  but 
not  so  high  as  to  degrade  the  probability  of  detection  below  a  tolerable 
level. 

The  last  two  factors  depend  almost  entirely  on  a  judicious  choice 
of  a  signal  waveform.   Hence  it  is  topical  to  consider  a  measure  or  cost 
function  by  which  various  waveforms  may  be  compared,   Siebert's-*-  modifi- 
cation  of  Woodward's   ambiguity  function  provides  such  a  comparison. 


11 


It  is  for  this  reason  that  the  ambiguity  function  forms  the  basis  of 
the  modern  radar  detection  philosophy. 

This  thesis  is  primarily  concerned  with  a  class  of  codes,  called 
complementary  codes,  which  minimizes  range  ambiguities .   Most  codes, 
which  have  been  described  and  analyzed  in  the  literature  on  radar  wave- 
form design,  have  the  drawback  that  range  and  doppler  ambiguities  exist. 
Examples  of  these  codes  are: 

1)  Linear  FM  (chirp) 

2)  Truncated  pseudo-noise  codes 

3)  Barker  codes. 

Of  the  codes  mentioned  above,  the  Barker  codes'  autocorrelation  functions 
exhibit  the  smallest  range  sidelobes,  and  have  been  named  "perfect"  words. 

Complementary  codes  have  the  feature  that  their  autocorrelation 
functions  have  no  range  sidelobes.   Although  much  work  has  been  done  in 
the  area  of  obtaining  properties  of  these  codes,  the  analysis  of  the 
ambiguity  functions  for  these  codes  has  not  been  done.   It  is  the  pur- 
pose of  this  thesis  to  analyze  complementary  codes,  in  both  range  and 
doppler,  for  two  cases: 

1)  Narrow-band  analysis 

2)  Wideband  analysis. 

By  narrow-band  we  mean  that  the  doppler  effect  does  not  significantly 
alter  the  envelope  of  the  transmitted  signal.   The  narrow-band  analysis 
will  be  formulated  by  utilizing  Siebert's  ambiguity  function.*   However, 


*Siebert*s  ambiguity  function  is  a  modification  of  Woodward's 
ambiguity  function.   The  essential  difference  is  that  Siebert's  ambiguity 
function  is  real  valued,  whereas  Woodward's  ambiguity  function  is  complex 
valued. 
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one  of  the  postulates  in  the  formulation  of  Siebert's  ambiguity  function 
is  that  doppler  docs  not  significantly  affect  the  envelope  of  the  trans- 
mitted signal.   Hence  Siebert's  ambiguity  function  cannot  be  used  in  the 
wideband  analysis.   We  will  extend  Siebert's  ambiguity  function  to  include 
the  wideband  case.   In  Appendix  B  we  will  define  the  region  of  validity 
of  the  narrow-band  analysis.   The  essence  of  the  derivation  of  the  region 
of  validity  of  the  narrow-band  analysis  is  that  the  time-bandwidth-velocity 
product  must  be  much  less  than  the  velocity  of  propagation  divided  by 
two,  in  order  that  we  may  assume  that  the  envelope  of  the  transmitted 
signal  is  not  altered  by  the  doppler  effect.   Then  by  wideband,  we  mean 
that  doppler  has  a  significant  effect  on  the  envelope  of  the  transmitted 
signal.   Throughout  this  thesis,  we  will  assume  that  target  acceleration 
is  negligible  over  the  interval  of  measurement. 

Let  us  now  consider  complementary  codes.   Complementary  codes  or 
series  were  originally  conceived  by  M.  J.  Golay,^*^  in  his  work  on 
infrared  multi-slit  spectrometry.   An  extensive  list  of  properties  for 
complementary  codes  was  formulated  by  M.  J.  Golay  and  additions,  based 
upon  applications  of  group  and  vector  theory,  were  formulated  by 
S.  Jauregui  in  his  doctoral  dissertation. 

A  set  of  complementary  series  may  be  defined  as  "a  pair  of  equally 
long,  finite  sequences  of  two  kinds  of  elements  which  has  the  property 
that  the  number  of  pairs  of  like  elements  with  any  one  given  separation 
in  one  series  is  equal  to  the  number  of  pairs  of  unlike  elements  with 
the  same  given  separation  in  the  other  series. ""  These  sequences  are 
binary  and  are  aperiodic. 

It  has  been  shown  that  the  autocorrelation  function  of  these  codes 
gives  rise  to  a  central  peak.   The  height  of  the  peak  is  dependent  on 
the  length  of  the  code,  for  zero  delay,  and  decreases  linearly  to  zero 
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for  delay  equal  to  the  duration  of  one  code  symbol,  then  remains  at  zero. 
It  is  this  property  of  the  complementary  codes  that  makes  this  coding 
method  attractive  for  use  in  the  detection  problem.   Because  the  defi- 
nition of  the  complementary  codes  appears  rather  formidable,  we  will 
utilize  the  property  given  above  whenever  we  refer  to  complementary 
codes.   Appendix  F  is  a  partial  list  of  the  properties  of  complementary 
codes.   This  list  is  provided  so  that  complementary  codes  of  almost  any 
length  may  be  generated  from  basic  length  codes  called  kernels.   A  kernel 
is  defined  as  a  basic  length  code  which  cannot  be  decomposed  into  shorter 
length  codes  by  an  inversion  of  the  standard  generating  methods.   To  the 
author's  best  knowledge  there  are  only  four  kernels  that  are  known  to 
exist.    These  kernels  are  listed  in  Appendix  F. 

At  this  point  xve  might  indicate  what  the  ideal  ambiguity  function 
would  be.   Since  we  would  like  to  resolve  two  targets  in  range  and  doppler 
to  an  infinite  degree,  the  ideal  ambiguity  function  would  be  a  delta 
function  at  the  origin  of  the  range,  doppler  axis.   Thusly,  the  syste- 
would  have  infinite  resolution.   This  can  never  be  achieved;  what  would 
be  required  is  a  code  of  infinite  length  with  the  duration  of  each  code 
symbol  inf initesimally  small.   However,  we  should  not  discard  the  idea 
of  an  ideal  ambiguity  function,  since  it  gives  us  a  reference  by  which 
to  compare  other  ambiguity  functions. 

A  brief  synopsis  of  each  chapter  will  be  now  given. 

In  Chapter  II  a  system  model  for  single  codes  is  postulated  and  a 
correlation  between  a  physical  waveform  in  the  system  and  the  ambiguity 
function  is  given.   By  single  codes  we  are  referring  to  binary  codes 
such  as  Barker  codes,  pseudo-noise  codes,  etc.,  i.e.,  codes  which  con- 
sist of  a  single  sequence;  as  opposed  to  complementary  codes  which 
consist  of  a  pair  of  equally  long  binary  sequences.   In  this  chapter 
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the  analysis  is  entirely  narrow-band.   The  single  code  ambiguity  function  is 
derived  to  ease  computational  matters  in  the  derivation  of  the  ambiguity  func- 
tion for  complementary  codes. 

Chapter  III  contains  the  derivation  and  the  justification  of  the  form  of 
the  narrow-band  complementary  code  ambiguity  function. 

In  Chapter  IV  a  system  model  is  postulated  for  the  complementary  codes 
and  the  derivation  of  the  complementary  code  ambiguity  function  is  accomplished. 
Computer  programs  were  run  for  the  two  kernels  of  length  10  and  the  kernel  of 
length  26,  and  ambiguity  diagrams  were  constructed.   One  property  of  the  entire 
class  of  binary  codes  in  which  phase-reversal  modulation  is  utilized,  is  un- 
covered in  this  chapter.   It  is  proved  that  in  the  doppler  domain,  the  zero 
crossings  of  the  ambiguity  function  move  closer  to  the  origin  of  the  range- 
doppler  axis  as  the  length  of  the  code  is  increased;  assuming  that  the  duration. 
T,  of  a  code  symbol  remains  constant.   The  implication  of  this  phenomenon  is 
that  better  doppler  discrimination  is  possible  for  longer  code  lengths,  all 
other  things  being  equal. 

Chapter  V  contains  the  investigation  of  the  doppler  effect  on  the  trans- 
mitted signal  and  the  proof  that  the  envelope  of  the  transmitted  signal  is 
altered  due  to  doppler.   This  chapter  in  conjunction  with  Appendices  A  and  D 
shows  that  the  doppler  effect  is  threefold: 

1)  The  carrier  frequency  is  shifted  by  the  doppler  frequency. 

2)  The  envelope  of  the  transmitted  signal  is  compressed  or  expanded 
according  to  whether  the  target  is  closing  or  opening  with  respect 
to  the  receiver. 

3)  The  received  and  transmitted  energies  are  different;  in  fact,  the 
received  and  transmitted  energies  differ  by  the  same  factor  as  the 
received  and  transmitted  carrier  frequencies.   The  above  statement 
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assumes  a  lossless  medium  and  ignores  spherical  spreading  of 
the  transmitted  and  received  waves. 

Chapter  VI  contains  the  derivation  of  the  wide-band  ambiguity 
function  for  single  codes  as  well  as  complementary  codes,  valid  for 
positive  differential  doppler  and  positive  range  deviation.   Unfortunately 
due  to  time  limitations,  computer  programs  to  generate  numerical  values 
for  the  wideband  complementary  code  ambiguity  function  were  not  accor- 
plished. 

To  clarify  the  discussion  of  the  model  of  the  system,  which  is 
postulated  in  Chapter  II,  a  brief  formulation  of  the  autocorrelation 
function  of  complementary  codes  will  be  given.   As  mentioned  previously, 
complementary  codes  consist  of  two  binary  sequences,  which  will  be 
referred  to  as  Code  A  and  Code  B  or  simply  A  and  B.   To  obtain  the  auto- 
correlation function'  for  the  complementary  code,  R  (7  ),  we  conpute 
the  autocorrelation  function  for  code  A,  RA(T),  and  for  code  B,  Rg(r). 
Then  RC(T)  is  obtained  by  summing  RA(7)  and  RB(T),  i.e., 

RC(T)  =  RA(T)  +  RB(T), 

The  property  of  interest  in  the  autocorrelation  function  of  the  comple- 
mentary codes  is  that  no  range  sidelobes  exist.   All  binary  codes  have 
autocorrelation  functions  with  range  sidelobes.   However,  RA ( T  )  has 
range  sidelobes  which  are  the  negative  of  the  sidelobes  of  RF(7). 
Hence  when  RA(T)  and  RB(T)  are  summed,  the  range  sidelobes  cancel.   Fig- 
ure 1-1   shows  this  unique  property  of  the  complementary  codes. 
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Figure    1-1 
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Autocorrelation  function  for  code  A,  length  10; 
A  =  set  J^-1,+1, +1,-1, +1,-1, +1,+1, +1,-1} 

Autocorrelation  function  for  code  B,  length  10; 
E  -  sfc-t  l-1  .+1,  +  1,M, +1,+1, +1,-1,-1, +1} 

Autocorrelation  function  for  the  complementary 
code,  l^r-;th  10. 
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CHAPTER  II 

MATCHED  FILTER  APPROACH  FOR  THE 
SINGLE  CODE  SYSTEM  (NARROW-BAND  ANALYSIS) 


It  might  be  well  to  look  at  the  single  code  system  before  attempting 
to  analyze  the  complementary  coded  system.   In  the  matched  filter 
approach  the  underlying  theme  is  to  correlate  the  ambiguity  function 
with  a  physical  waveform  in  the  system.   The  ambiguity  function  is  an 
excellent  measure  of  range  and  velocity  resolvability ,  i.e.,  how  closely 
one  can  distinguish  two  targets  in  range  and  velocity.   A  cornon  ris- 
conception  is  to  state  that  the  output  of  the  matched  filter  is  identical 
to  the  ambiguity  function  for  all  cases.   This  statement  holds  some 
truth;  however,  it  is  not  generally  correct  in  carrier  modulated  systems. 
Let  us  depart  from  this  discussion  to  model  our  system  before  attempting 

to  clarify  the  above  statement.   The  model  that  is  chosen  is  a  slight 

8 
modification  of  the  model  used  by  Davenport. 


Transmitter  j s(t)  J  Transmitted  \       s1 (t-td)  +  n(t)  J  Receiver 

'  MpHium 


Medium 
information:         'disturbance 


source 
(target) 


(additive) 


System  model 
Figure  2-1 

The  difference  between  Davenport's  model  and  the  model  that  is  chosen 
is  that  in  the  radar,  as  well  as  the  sonar  case,  the  information  source 
interacts  in  the  medium,  whereas  in  Davenport's  system  model,  the  infor- 
mation source  interacts  in  the  transmitter. 
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The  transmitted  signal  is  chosen  as  a  phase-reversal  modulated 

signal,  where 

s(t)  =   \cos(u  0t  -  n.p   )  2-1 

where  n^  =  j  0    depending  on  the  value  of  the  ith  code 

'  1 

symbol 

\  =  normalization  constant . 

For  convenience  the  signal  energy  is  normalized  to  unity  by  specifying 
that 

ls2(t)  dt  =  1 

It  is  a  well  known  fact  that  a  phase-reversal  modulated  signal  is 
identical  to  a  DSB/SC  signal;  hence  the  transmitted  signal  can  be  written 
as 

s(t)  =  I(t)  cos  u0t,  2-2 

where  I(t)   ■  V    ~>.    \i  -      (t-iT)   -    Jx  -    (i  +  1)TJ\ 

+       d 
and  X£  =  -  1,  depending  on  the  value  of  the  ith  code 

symbol. 

The  normalization  constant  can  be  easily  evaluated  by  squaring  equation 
2-2  and  performing  the  integration.   Thus 

,2    NT 
s2(t)  dt  =  -j-   \  \l  +   cos2,  0t!  dt  2-3a 

2  ? 

=     x   NT  cXzBin2t0oNT  2-3b 

2  4ti0 

The  second  term  in  equation  2-3b  is  very  small  compared  to  the  first 
term,  for  L0_  large.   This  term  can  be  interpreted  as  the  difference 


19 


between  the  initial  and  final  value  of  the  carrier  phase  and  trill  he 
snail.   Hence  the  normalization  constant  is  approximately  given  by 


o< 


2  o  , 

— —  I  —  3c 

NT  « 


The  received  signal  contains  a  noise  term,  n(t),  which  will  be 

9 
ignored  since  the  receiver  will  utilize  a  matched  filter   to  raxir.ize 

the  signal-to-noise  ratio.   Thus  the  received  signal,  s'(t-t,),  is  the 
transmitted  signal  modified  by  the  information  source,  and  contains 
range  and  velocity  information.   Let  us  assume  that  the  target  speeds 
are  sufficiently  less  than  the  speed  of  propagation  of  the  wave  in  the 
medium  (say  of  the  order  of  10   x  c,  where  c  is  the  velocity  of  propa- 
gation,)  The  principal  doppler  effect,  in  this  case,  is  to  shift  the 
carrier  frequency  by  the  amount  f  ,,  where  f ,  is  the  doppler  frequency; 
doppler  will  have  little  effect  on  the  duration  of  a  code  symbol,  T. 
(See  Chapter  IV  for  the  discussion  of  the  doppler  effect  on  the  duration 
of  a  code  symbol.)   With  the  above  conditions  postulated,  the  received 
signal  may  be  given  by 

s'C-fc-tj)  =  Kt-V>  cosftuvvu^l  v 

=    °<2X.{i(-t,-ta-5T)-l[t-td-(L  +  i)Tlco^; .  .     ftt-ty],         2-4 

where  the  convention  is  chosen  so  that  the  doppler  frequencv  is  positive 
for  closing  targets,  and  negative  for  opening  targets. 

At  this  point  the  question  might  be  raised  as  to  what  are  the 
various  waveforms  in  the  receiver.   As  mentioned  earlier  in  this  chapter, 
a  matched  filter  will  be  utilized  in  the  svstem  model.   To  obtain  the 
signals  to  correlate  with  the  received  signal,  we  translate  the  trans- 
mitted signal  by  various  expected  doppler  frequencies  as  indicated  in 
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Frequency    shirter 


si  (t) 
T(t-tri)  „ 


■  i(t) 


sn(t) 


f— — i * Xi^mo  > 


Dorplf  r 
fi Iter 

bank 


'  XCr,^)  » 


'•  X(rt#A«>i 


Comparator 


to  display  units- 


Figure  2-2a 


r, 


s'( 


t-td) 


i(t) 


l^r — 
•E — ( 


£ — »|Integrator 


M<r,)W)a;) 


Collator  - — i 


a  a 


Lm 


-^Integrator!- 


^Crv^t^ao 


Integrator 


4(rW)u3ai> 


tc  Comparator  •* 


X(r,  u>ac) 
4 


"Envelope 
detector" 


Figure    2-2b 

Figure  2-2a   Receiver  block  diagram 

2-2b   ith  doppler  filter  of  the  filter  bank 
in  Figure  2-2a 
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Figure  2-2a.   The  doppler  shifted  signal  is  then  delayed.   The  delay  is 
required  to  account  for  the  propagation  time  to  and  from  the  target.   To 
instrument  this  system  in  real  time,  the  delay  mechanism  might  be  of  the 
form  of  a  tapped  delay  line.   We  shall  denote  the  delayed,  doppler  shifted 
signal  as  the  replica  signal.   The  replica  signal  is  given  by 

si(t  "  *k)  -  l'(t  -  Tk)cos  [(i0o  +  (Adi)(t  -  rk)]  2-5 

Let  us  now  analyze  the  various  signals  in  the  receiver.   To  simplify 
the  analysis,  the  matched  filter  block  diagram  is  duplicated  from 
Figure  2-2b. 

The  output  of  the  matched  filter  is 

yk(rk^'>"W>  =  \  S'a-ta)S;tt-^<H  2-6a 

-co 


Ki-r^ ccs^lo^u-t^  <H  2"6b 


Then  let 


fcl  =  t  -  td  2-7a 

&ox    =   ^o  +  ^di  2-7b 


Equation  2-6b  becomes 


,co 


Cos  [  ^(-t^-la  -  vl  <&!   .         2"8 
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Then  let 

X     -rk-td  2.9a 

wd  =  u>d  -  o)di  .  2-9b 

If  the  above  changes  in  variables  are  effected,  then 


oo 

yk(T,u>d)   =   (ia^Kvt)  cosC^  +  tAiaHi  cos  u)0i(vO<ft. 


2-10 


The  variables  X   and  ££),  deserve  some  comment.   From  equation  2-9a, 

we  see  that  if  X    is  zero,  then  the  propagation  delay  to  the  target 

Tk 
has  been  matched  in  range  bin  T,  .   Hence   —  is  the  range  to  the 

k  2 

target.   From  equation  2-9b  we  see  that  if  06d  is  zero,  then  U3d-^  is 
the  doppler  frequency.   Thus  both  X.    and  u)j  represent  the  departure 
of  the  target's  range  and  velocity  from  the  estimated  range  and  velocity. 
The  variable  t«  is  just  the  variable  of  integration,  and  since  it  is 
inconsequential  as  to  what  nomenclature  is  ascribed  to  it,  we  shall 
denote  t^  by  t.   00oi  represents  the  estimated  carrier  frequency  of  the 
replica  signal.   Hence  the  output  of  the  matched  filter  in  the  syster 
is  given  by 


00 

r 


%(.?,^     =      \  KV)I(t-r)   cos  (  i001+ uidH   cos  u)0_L(-t-T)cit    ^  2-lla 

-oo 
=    z\Itt)I(^T)    COS((jOdt-LG01r)dt 

-Oci 

■+l\lCt)IC-t-r)cos[(2w0i-»-uiid)t-ui01t:ldt   .  2-iib 
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The  second  term  in  equation  2-llb  has  frequencies  centered  about  twice 
the  carrier  frequency  plus  the  doppler  frequency,  and  will  be  rejected. 
Note  that  the  output  of  each  matched  filter  is  an  ordinate  in  the  *C 
domain. 

Next  we  require  a  device  to  collect  the  output  of  each  of  the 
matched  filters  in  an  orderly  fashion.   Hence  the  use  of  the  term, 
collator,  to  describe  that  device  is  appropriate.   We  shall  define  the 
collator  in  Figure  2-2b  as  that  device  which  collects  the  various  out- 
puts of  the  matched  filter  sequentially  as  a  function  of  range  deviation, 

.   To  clarify  the  above  definition,  we  might  think  of  the  collator  as 
a  graph  plotter,  where  ordinates  (outputs  of  the  matched  filters)  are 
plotted  as  a  function  of  range  deviation,  T  .   The  "envelope  detector" 
in  Figure  2-2b  is  not  an  envelope  detector  in  the  normal  sense,  and 
requires  some  explanation.   As  mentioned  above,  the  collator  collects 
the  outputs  of  each  matched  filter  and  produces  a  continuous  curve  in 
the  X.   domain.   This  curve  oscillates  at  i0ol  in  the  X   domain,  but  only 
its  envelope  has  information  content.   The  "envelope  detector"  refers 
to  that  device  that  takes  the  envelope  of  the  output  of  the  collator, 
but  in  the  range  deviation  or  T  domain.   In  the  usual  sense  of  the 
term,  envelope  detection  is  accomplished  in  the  time  domain  and  hence 
a  distinction  must  be  made  in  our  usage  of  the  term  envelope  detector. 

Now  let  us  look  at  the  output  of  the  collator; 

CO 

UlT.lBil  =  A{    Kt)l(t-r)  o-;(co.t-w.,r)dt  .  2-12 


If  the  Euler  identity  for  the  cosine  function  is  utilized,  the  output 
from  the  collator  is 
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LJV  \   _..-.,.    ^A^ 


-oo 

^i^'lX\  lU>IU-r)eil"jtdt    .  2-13 


4 

-oo 


If  we  note  that  the  second  term  of  equation  2-13  is  the  complex  conjugate 
of  the  first  term  and  that  the  sum  of  a  quantity,  say  z,  and  its  corplex 
conjugate  is 

Z  +  Z*  =  2  Re(Z)  , 

then 

CO 

jjCe,^)-  Z(J?e{iei^\lWia-r)e5^at}.       2-14 

-do 

The  terms  I(t)  and  I(t  -  T  )  are  the  modulation  functions  and  involve 
the  limits  of  integration  and  the  values  of  the  codes.   The  modulation 
functions,  according  to  equation  2-2  are 


I0U  =  ^ZX;fl^-iT)-Mt-lM)T"  2-15a 

1=1.    *-  J    * 

M 

±(i-o  =  o(2xt(fli4-T:-if)-llU-,c-e+iVtli(  2-i5b 


The  unit  step  functions  in  equations  2-15a  and  2-15b  will  determine  the 
limits  of  integration  and  the  Xi's  are  ±1,  depending  on  the  value  of  the 
codes.   Figure  2-4  is  the  plot  of  the  time  waveforms  of  I(t)  and  I(t  -  r  ) 
To  obtain  the  limits  of  integration,  we  merely  observe  the  zero  crossings 
of  the  functions.   Hence  the  output  of  the  collator  is  given  by 
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The  above  equation  can  be  generalized  to  a  code  of  any  length  Na  and 
will  be  written  as  the  sum  of  two  series.   Therefore  the  output  of  the 
collator  is 


M 


^,^  =  f  ^(e^g^Xk,!^    e»"*\H 


i=  1  J 

CKm)T 


where  N  =  length  of  the  code 

k  =  smallest  integer  >  \  T/t  \ 

T  =  the  duration  of  a  code  symbol  =  clock  period 

The  first  integral  in  equation  2-16  can  be  evaluated  as 


CA  =    \     e3     at    = z_§ 

Tt(i-1)T  3  d 


2-17a 
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=  -j^  -[cos  u^Ck+i)T  +  j  sm  u^(H+iVT  -  co6tod[r+(i-i)T] 
-Jsinu^Lt  +  Ci-DT]  "^  e 

If   the   following   trigonometric   substitutions  are  utilized, 
cos  x-  cos  y    =    -  Z  sin  (*^-)  sm  (^-) 

equation  2-17b  becomes 


C,  =  ^  [-z  sm  S>Pa*fc-"T*<i  sm  ^^-r] 


i    .uur(2uk-Di->g]  -, 


2-17b 


2-18a 


2-l8b 


The  second  integral  in  equation  2-16  can  be  evaluated  as 


G2=  ^   e    di  -  —    ,^  .  2-19a 


By  similar  arguments  used  to  derive  equation  2-18d,  C„  becomes 


C^i-Sia^Me      *         .  2-1% 
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We  then  substitute  equations  2-18b  and  2-19b  into  equation  2-16. 

e**"^:  ^^  2-20 


If  the  exponential  terms  are  combined*  and  the  real  part  of  equation 
2-20  taken,  then 


m  p  ^Cr-kT)  c^  r^(2uk)T4C2^+iAQt:l    2_n 


2 


In  order  to  plot  equation  2-21,  it  is  necessary  to  obtain  the  equation 
in  the  form  of  an  envelope  term  multiplied  by  an  RF  terra.   Let 


.  ^(auk-i)T 

% 7T ^  2-22a 


*,-    frC21+k)T  2.22b 

2*  ^ 


•00,  =     2<A»H-<4i  2-22c 

1         z —      > 
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E2=  XiXk^  s,„  i^ML       # 


2-22d 


2-22e 


Then 


2-23a 


We  now  utilized  the  trigonometric  identity 


cos (a  +  b)  =  cos  a  cos  b  -  sin  a  sin  b  , 


Equation  2-23a  becomes 


ifCc,^)  =  —  j2  Ej.  cos  a,,  cos  aAt:  -  2  E^in^siniitf: 


ikk+i) 


~^U<Z  E.eosi,  -^  E^gosxJ  cos  li^c 
2  EAsin^  +A  E^siniJ  sin  u\r\  . 


2-2  3b 


Then 


A  cos  a  -  B  sm  a   = 


A 


A2  +  B2      cos  (a-  -W-|-) 
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Thus   the   output   of    the   collator  becomes 


est1  frK"k  ^S  -r 


ri-k 


^I-(k4i) 


}k 


+  [21  Ej.siii^  +  2  E,  sin  ^f\    c°5teaT-0) 


2-23c 


where  0  is  the  phase  due  to  the  amplitude  of  the  sine  and  cosine  terms. 


If  equations  2-22  a,  b,  c,  d,  e,  are  substituted  into  equation 

2-23c,  and  if  the  sine  terms  are  placed  in  the  familiar  sin  x  form, 

x 

then 
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Equation  2-24  is  plotted  in  Figure  2-5;  note  the  RF  fine  structure 
under  the  envelope. 

Now  let  us  analyze  the  ambiguity  function.   Essentially,  definitions 
and  results  will  be  stated  from  Siebert's  paper  with  the  exception  of 
a  few  terms  which  will  be  modified  to  conform  to  the  nomenclature  of 
this  thesis.   The  ambiguity  function  as  defined  by  Siebert  is 


4Cr,u^-^r 


&(t)S0*tt-r)e1     & 
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where  S  (t)  is  the  complex  envelope  of  the  transmitted  signal. 


In  terms  of  this  thesis 


S(t)  =  I(t), 


I(t)  is  real.   Therefore 


*(«.<*)  =i 
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Equation  2-26  departs  from  equation  2-14  in  that  the  carrier  phase  term 
is  absent  and  the  magnitude  of  the  integral  is  taken.   Utilizing  the 
same  analysis  used  in  obtaining  the  limits  of  integration  for  the  output 
of  the  collator,  we  see  that  equation  2-26  can  be  written  as 
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If  the  results  of  equations  2-17  through  2-20  are  utilized,  equation 
2-27  becomes 


4(*,«tf  =  &■ 
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3  ^  (2Uk-DT  5  ^ 


N-M 


^    "e  -  -v  ^xp^ 


2> 
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m 


The  term   e  a"   can  be  factored  from  equation  2-28.   Note  that 


e1  ? 


-   i 


If  the  magnitude  of  equation  2-28  is  taken,  then  the  enuation  in  its 
final  form  is 
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— -2 — 

2  .  »/Z 

sin  u}a(2uk)>"  I-  \  2-29 


In  the  earlier  portions  of  this  chapter,  the  statement  was  made 
that  the  output  of  the  matched  filter  is  identical  to  the  ambiguity 
function.   We  are  now  prepared  to  discuss  this  particular  point.   In 
comparing  equation  2-24  (the  output  of  the  collator)  with  equation  2-29 
(the  ambiguity  function),  we  note  that  the  ambiguity  function  is  identical 
to  the  envelope  of  the  output  of  the  collator,   hence  the  difference 
between  the  ambiguity  function  and  the  output  of  the  matched  filter  is 
that  there  exists  an  RF  fine  structure  of  frequency  f   ,  in  the  output 
of  the  matched  filter  which  is  absent  in  the  ambiguity  function.   Note 
that  if  the  phase  of  the  RF  fine  structure  is  some  odd  multiple  of  -*-, 
then  the  output  of  the  matched  filter  will  be  zero.   Thus,  we  can  say 
qualitatively  that  the  ambiguity  function  measures  the  positive  raxirurs 
of  the  output  of  the  matched  filter  per  cycle  of  the  RF  fine  structure. 
To  circumvent  the  problem  of  outputs  of  the  matched  filter  being  zero 
near  the  central  peak,  we  sum  and  take  the  envelope  of  the  output  of 
the  matched  filters.   However,  care  must  be  taken  in  choosing  the  differ- 
ential range  deviations.   Suppose  that  the  differential  range  deviation, 
£X     ,  is  some  multiple  of  the  inverse  of  twice  the  frequency  of  the  RF 
fine  structure,  i.e., 
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Furthermore,  suppose  that  in  any  matched  filter,  the  phase  difference 

between  the  transmitted  signal  and  the  received  signal  is  an  odd  multiple 

of  -3L-  radians.   Then  the  output  from  all  the  matched  filters  will  be 
2 

zero.   This  point  is  easy  to  see  because  the  system  will  be  tracking  the 
zero  crossings  of  the  RF  fine  structure. 

The  RF  fine  structure  is  also  troublesome  in  the  area  of  multi- 
plexing.  This  point  will  be  discussed  in  Chapter  IV,  when  we  will 
formulate  ideas  on  the  type  of  multiplexing  to  be  used  in  the  complemen- 
tary coded  system. 

The  output  of  the  "envelope  detector"  is  given  by  the  envelope  of 
equation  2-24.   However,  as  mentioned  earlier,  the  output  of  the  "enve- 
lope detector"  is  identical  to  the  ambiguity  function. 
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CHAPTER  III 

DEFINITION  OF  THE  AMBIGUITY  FUNCTION 
FOR  COMPLEMENTARY  CODES 


In  this  short  chapter,  we  shall  extend  Siebert's  Ambiguity  function 
to  the  ambiguity  function  for  complementary  codes.   Let  us  analyze  the 
signal  processing  techniques  required  to  produce  the  desired  results  in 
the  complementary  coded  system.   In  Chapter  II,  it  was  shown  that  in 
the  single  code  system,  the  output  of  the  "envelope  detector"  in  Figure 
2-2b  was  identical  to  the  ambiguity  function.   What  then,  should  be  the 
configuration  of  the  complementary  coded  system?   Should  the  complementary 
coded  system  be,  in  essence,  two  single  coded  systems  such  that  the 
outputs  of  the  two  "envelope  detectors"  be  summed  together?   Or  should 
the  two  codes  be  summed  elsewhere? 

Let  us  postulate  the  complementary  coded  systerr  to  be  comprised  of 
two  identical  single  code  systems  with  the  outputs  of  the  "envelope 
dectors"  of  code  A  and  B  summed  together.   In  this  case  the  ambiguity 
function  of  the  complementary  coded  system  will  be  the  sum  of  the  ambi- 
guity functions  of  code  A  and  B,  i.e., 


±zLr^)  =  %(x^6)  +!i BCr,utf  . 


3-1 


One  of  the  properties  of  the  ambiguity  function  is  that 


4(r^ 


-  0( 
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Rix) 


Hence  equation  3-1  for  \x>  ^  =  0  becomes 

As  indicated  in  the  introduction  of  this  thesis,  the  desired  property 

of  the  complementary  codes  was  that  off  the  central  peak,  the  autocorrelation 

function  of  code  A,  R^CT),  was  the  negative  of  the  autocorrelation  function 

of  code  B,  RgCf),  so  that  when  the  autocorrelation  functions  were  summed, 

the  range  sidelobes  canceled.   Since  the  proposed  definition  (equation 

3-1)  of  the  ambiguity  function  contains  two  positive  quantities  |R^(t)\   , 

and  JRb(T  )|   »  the  definition  is  inadequate. 

As  a  second  approach  to  the  problem,  let  us  define  the  ambiguity 
function  of  the  complementary  coded  system  as 


^°  \*& 


*«.(** <Ai  =  \ \\  IfcMljjfc-rte1    '  &> 


* 


3-3 


where  I^Ct)  ~  envelope  of  code  A 
Ig(t)  =  envelope  of  code  B, 

For  zero  doppler  frequency,  equation  3-3  becomes 


±t(X)Q)m   OlA&A0c)+R*M\.  3-4 
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It  is  easy  to  see  that  the  range  sidelohes  of  code  A  will  cancel  the 
range  sidelohes  of  code  B  so  that  the  range  sidelohes  for  the  connle- 
mentary  coded  amhiguity  function  will  he  nonexistent  for  zero  doppler 
frequency.   Equation  3-3  will  he  chosen  as  the  definition  of  the  ambig- 
uity function  for  the  complementary  codes. 

Equation  3-3  gives  us  the  clue  as  to  the  configuration  of  the 
complementary  coded  system.   Ue  can  infer  that  the  two  codes  should  be 
summed  at  the  output  of  the  matched  filters,  then  envelope  detected. 
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CHAPTER  IV 

MATCHED  FILTER  APPROACH  FOR  THE  COMPLEMENTARY 
CODED  SYSTEM  (NARROW- BAND  ANALYSIS) 


A .  BACKGROUND 

In  the  complementary  coded  system,  some  form  of  multiplexing  is 
required  since  two  codes  will  be  utilized  to  obtain  the  desired  cancel- 
lation of  the  range  sidelobes.   Of  the  three  obvious  types  of  multi- 
plexing (time  division,  frequency  division,  and  quaternary  coding)  only 
time  division,  TDM,  and  frequency  division  multiplexing,  FDM,  v;ill  be 
discussed. 

B.  PROPOSED  FDM  SYSTEM 

Assume  for  the  moment  that  the  FDM  system  is  modelled  as  below. 
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Complementary  coded  FDM 
transmitter  block  diagram 

Figure  4-1 


It  might  be  well  to  observe  the  transmitted  spectrum  to  aid  in  the 

understanding  of  the  signal  processing  schemes  that  will  be  required  in  the 
receiver.   The  frequencies,  f-^  and  f~,  are  at  RF.   The  outputs  of  each 
of  the  phase-reversal  modulators  can  be  easily  obtained  from  equation  2-2 
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if  appropriate  subscripts  are  used.   Hence 


1=1 


>Buw  oic2XBf(lltt-irt-llit-(i4ijr]  ccs^  . 


i*i 


4-2 


It  is  easy  to  see  that  the  spectrum  of  the  spectrum  of  the  transmitted 


signal  consists  of  two 


sin  x 


functions  centered  at  f-^  and  f2 


One-sided  spectrum  of  the  transmitted  signal 
Figure  4-2 

Thus  to  detect  each  code,  the  system  requires  bandpass  filters  centered 
at  the  carrier  frequency  plus  the  doppler  frequency.   Fence  the  system 
block  diagram  is  given  by  Figures  4-3a  and  4-3b.   Since  the  output  of 
each  range  bank  can  be  treated  the  same  as  in  the  single  code  syster, 
equation  2-15  applies  with  the  substitution  of  appropriate  subscripts. 
Let  &  be  the  envelope  quantity  in  equation  2-15;  hence  the  output  from 
any  range  bank  in  Figure  4-3b  for  code  A  is  given  by 


42 


E 

03 

U 

hi) 

ca 

•H 

T3 

CTS 

r~\ 

J4 

l 

V 

-J- 

O 

r~\ 

CO 

X 

U 

3 

U 

tun 

0) 

•H 

> 

&H 

•H 

CU 

o 

<u 

fc 

s; 

Q 

fe 

O 
P 

03 

Jh 

03 

E 
O 

o 


CU  ► 
a  Jh 

o  o 

rH  p 

CU  o 

>  0) 

C  -p 

U3  d) 

t  t3 


3 
CO 


CO 


•3 


o  u 

C    CU 

CU    -P 

ct-h 
cu  xi 


CO  O 


e* 


CD 

C  C 
03  03 
CtJ  X) 


rr 


\.s- 


0) 

p 


CQ 


«H 


V 

CO 

p 

•H 
C 

>> 

03 
rH 
CX 
co 

•H 

O 
P 


0)  c 

o  o 

-H  p 

CU  o 

>  0) 

C  P 

W  CD 

c  T3 


3 

CO 


an 


CU 

c  c 

03    03 
CtJ  JD 


TT 


Oh 

pq 


3 


CD 


CU 

c  c 

03     (D 


CO 


CU 
cc   cO 

ct;  xi 


IT 


3 


""S^ — ;=: 


C-, 

CU 

p 

PQ  -h 


CU  P 
3  Vh    0) 

CD  x:   o 
^  co  o 


T3 

p 

I 
p 


co 


43 


3 
10 

o 
p 


2 


o 
-p 


O 

o 


4D 

<^\ 

I 

CD 

W 


03 

I 

hO 
•H 
fin 

c  < 


0) 

o 


03  O 

CD 
hQ 
C 
03 

x: 
p 

•H 


< 
to 


4 

7T* 


4. 


o 
p 

03 

u 

fafi 
0> 
P 

c 


'_ 
c 
p 
td 
t- 

o; 
p 

C 


»fW 


G        <u 

O  PL,    p 

t-t  CQ  rH 

Cm        -H 

Ch 


A  A 


The  output  from  the  corresponding  range  hank  for  code  B  is 

The  output  from  the  summer  in  Figure  4-3a  is  the  sum  of  the  outputs 
from  the  range  hanks  of  code  A  and  B.   Therefore 

The  effect  of  the  different  frequencies  of  the  RF  fine  structure  for 
code  A  and  B  is  to  lower  the  output  from  the  matched  filter  at  the 
central  peak  and  to  increase  the  range  sidelobes.   It  is  true  that,  if 
the  carrier  frequencies  of  both  codes  are  nearly  equal,  the  range  side- 
lobes  will  be  very  snail;  and  in  fact,  might  be  quite  tolerable.   However, 
bear  in  mind  that  the  characteristic  of  interest  in  a  complementary 
coded  system  is  for  range  sidelobes  to  be  nonexistent  for  zero  doppler. 

In  this  system  the  factor  which  is  troublesome  is  the  RF  fine 
structure  which  is  absent  in  the  ambiguity  function.   Hence  the  analysis 
provides  a  further  contradiction  to  the  statement  that  the  output  of  the 
matched  filter  is  identical  to  the  ambiguity  function.   Indeed,  it  is 
that  factor  which  is  different  between  the  output  of  the  matched  filter 
and  the  ambiguity  function  which  defeats  this  method  of  multiplexing. 

45 


C.   TDM  SYSTEM 

Another  multiplexing  method  is  TDM.   The  suceeding  portions  of 
this  thesis  will  dwell  entirely  on  phase-reversal  TDM  systems  for  comple- 
mentary codes.   The  results  from  Chapter  II  concerning  single  code 
systems  can  be  easily  extended  to  the  complementary  coded  systers.   In 
the  transmitter  the  essential  difference  is  the  time  multiplexing  circuitry; 
hence  the  transmitter  can  be  modelled  as  below. 
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TDM  transmitter  block  diagram 
Figure  4-4 

As  in  Chapter  II  the  transmitted  energy  will  be  normalized  to  unitv 
by  specifying  that 


DO 

^  sz(V)  <fe  .  1  . 


4-5 


To  obtain  the  normalization  constant,  we  note  that 


S(t)    = 


SA("t)     y        ior  0  4*<   HT 


CHW    , 


HT<^<  ZHT 
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sA(t)   and   sB(t)   are  given  by 


sAW  =   C^2xAi{A(VST)-Jl[t-(i+i)T]c1osu5o4      ,  /^6a 


;here  XA •  =  _1,  depending  on  the  value  of  the  ith  code  symbol  of  A 
XB^  =  ll,  depending  on  the  value  of  the  ith  code  symbol  of  B. 


Then 


qo  f>  9° 


\4A&(a  =   o(ciiI  +  ac  sul5&SI  «  otm 
z  4^  a 


>      B  ?  o   .  9 


Therefore   the  normalization  constant,      o^c»    -"-s  aPP*"oximately 


ac  = 


A 


HT 


4-7 


In  the  transmitter,  code  B  is  delayed  and  summed  with  code  A,   The 
receiver  is  modelled  in  Figure  4-5.   To  detect  the  signal  the  gating 
circuit  in  the  receiver  separates  the  two  codes.   The  output  from  each 
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TDM  receiver  block  diagram 
Figure  4-5 


48 


correlator  in  Figure  4-5  is  identical  to  the  output  from  a  correlator 

in  the  single  code  system.   Hence  the  output  from  the  "envelope  defector" 

is  given  hy 


X(r,U3d)  =  Envelope  of  {y^^ -v  y^T,^)}  , 


4-8 


Upon  examination  of  equation  2-15,  we  find  that  the  only  terms  which 
depend  upon  the  value  of  the  codes  are  the  X-'s.   The  rest  of  the  terms 
are  common  to  both  yA  (  T  ,  L&d)  anc*  Yb  (  *C »  ^d^'   hence  if  the  signed 
products  in  equation  2-15  are  replaced  by  signed  products  for  code  A  and 
B,  i.e. , 

X^X-  is  replaced  by  (XA^XAj  +  XB-jXBj) 

and  if  the  envelope  of  equation  2-15  is  taken >  then  the  output  of  the 
"envelope  detector"  will  be  given  by 


Z 
-rXBjXe^:]  coau^UUk-m-  +  (t-kr)  Sin  ^Jd  ^ 
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As  stated  in  Chapter  II  the  output  from  the  "envelope  detector"  is  iden- 
tical to  the  ambiguity  function.   Therefore  the  ambiguity  function  for 
the  complementary  codes  will  be  given  by  equation  4-9. 

Let  us  now  analyze  the  narrow-band  ambiguity  function  for  comple- 
mentary codes.   For  zero  doppler  equation  4-9  reduces  to 


+ 


(r-i. 


2 


1=1 


=     afji^to  +J^(-o|  ( 


4-10 


Thus  we  see  that  from  the  analysis  done  in  the  introduction,  the  ambiguity 
function  for  zero  doppler  contains  no  range  sidelobes.   We  note  that 
this  property  is  a  direct  result  of  the  RF  phase  synchronization  assumed 
in  equations  4-6a,  b. 
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For  the  case  of  zero  range  delay,  the  ambiguity  function  reduces  to 


n^iZ  * 


2v 


-n^  &*-£W, 


+  X6i.X5ila,^LlK^ 


/. 


4-U 


Note  that 


X.X.  =1  . 
1  l 


Then  if  the  sine  terms  are  factored,  equation  4-11  reduces  to 


4-12 


If  the  Euler  identity  is  utilized 


e  =  cos-e-  •+  jsm-©- 


and  noting  that 


i/ 


then,  we  see  that 


|e*^|  =  le^e^'T|=||e^' 


2,  Cos  uV4T   "^  \  2  ^in  0VtT 
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Thus  equation  4-12  can  be  written  as 


%(o^  *   Ok 


2. 


i  2  cos  ^T  +  i  J  siji  ,Vri 


4-13 


10 


If  the  following  trigonometric  formulas^  are  utilized 


£  cos  10  -  cos^^  cosec^ 


f  sinMy  ^   ^i&fi^^r.^ 


1=1 


a   C05ec'* 


and 


2_\  cosiS-  +  h  Za  sin  10  = 


thus  equation  4-13  becomes 


3^0,^=  ao 


T81ni¥ 


2. 


jjMn.%iI  oos  !Aff  es     * 


4-14 


Then  noting  that  the  absolute  value  of  a  product  is  the  product  of  their 
absolute  values,  and 

Finally  we  see  that  equation  4-14  becomes 


5^(0,^=  < 


KTsln 


2x 


2, 


4-15 
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If  equation  4-15  is  plotted,  one  interesting  result  is  observed 
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Ambiguity  function  for  complementary  codes} 
in  the  doppler  domain. 
Figure  4-6 

The  first  zero  crossing  of  the  ambiguity  function  occurs  at  «i_.;  succeeding 

1 
zero  crossings  occur  at  multiples  of  — -.   Kence  the  coded  pulse  sequence 

NT 

is  acting  like  a  single  pulse  of  length  NT  in  the  doppler  domain.   For 
a  better  insight  into  this  phenomena,  let  us  review  the  origin  of  the 
sin  •  '  "'   term.   A  glance  at  equation  4-12a  shows  us  that  the  term  origi- 
nated as  a  sum  of  exponentials,  i.e., 


if   .Utt(2L-i)T 
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1=1 
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The  criterion  for   the  ambiguity   function  to  be   zero   in  the  doppler  domain 
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Or  alternatively, 


-tiu^T  =    ii2.iT  •,       ja  =  A,2,  ■  ■■ 


4-17b 


Hence  we  see  from  equation  4-16  that  when  the  total  phase  shift  is  a 
multiple  of  2T  ,  the  ambiguity  function  for  zero  range  delay  will  be 
zero.   This  effect  can  be  viewed  as  the  addition  of  phasors  such  that 
when  the  function  is  zero,  the  phasors  have  added  to  forir  a  closed  loop. 
This,  of  course,  occurs  when  the  total  phase  shift  is  some  multiple  of 
2tt  . 

One  added  comment  should  be  made  concerning  the  zero  crossings  of 
the  ambiguity  function  in  the  doppler  domain.   Since  the  zero  crossings 
occur  at  iff,    then,  as  the  length  of  the  code  increases,  the  ambiguity 
function  will  be  compressed  closer  to  the  origin;  hence  approaching  the 
ideal.   However,  the  length  of  the  code  cannot  be  increased  indefinitely 
without  complications.   It  will  be  shown  in  Chapter  VI  that  longer  codes 
will  be  more  sensitive  to  changes  caused  by  the  doppler  effect,  and 
unfortunately  the  length  of  the  code  cannot  be  increased  indefinitely 
without  reducing  the  output  from  the  matched  filter.   The  above  compression 
effect  on  the  ambiguity  function  is  not  a  property  possessed  solely  by 
the  complementary  codes.   For  example,  single  phase-reversal  codes  have 
ambiguity  functions  which  experience  the  compression  effect  for  long 
codes.   It  is  the  length  of  these  codes,  rather  than  their  internal 
structure,  which  produces  compression  in  the  doppler  domain;  assuming 
that  the  duration  of  a  code  symbol,  T,  is  held  constant. 
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D.   SUMMARY  OF  RESULTS 

Equation  4-9  was  coded  into  Fortran  and  programs  were  run  for  the 
two  kernels*  of  length  10  and  the  kernel  of  length  26.   Plots  of  the 
ambiguity  function  were  made  for  constant  "doppler  phase./  $,  »  where 
this  is  defined  as  the  recurrent  argument  in  our  equations  for  the  ambi- 
guity function,  i.e., 


A 


where  T  =  duration  of  a  code  symbol. 

To  convert  doppler  phase  into  radial  velocity,  we  recall  that  for 
the  narrow-band  case 


Therefore, 


T  _  C  "£<!   _   c   §k  .  ... 

v_  =  -7T  ~r-      -  ~5 — ; 4-19 


Table  1  gives  the  conversion  from  doppler  phase  to  radial  velocity  for 
various  values  of  carrier  frequency  and  duration  of  a  code  symbol. 
Doppler  phase,  ^i  ,  was  defined  in  equation  4-18^  so  that  the  ambiguity 
function  need  not  be  recomputed  for  different  values  of  f  and  T. 


*As  stated  in  the  introduction,  a  kernel  is  a  basic  length  code 
which  cannot  be  decomposed  into  shorter  length  codes  by  an  inversion  of 
the  standard  generating  methods.   Conversely,  complementary  codes  of 
longer  lengths  can  be  generated  by  appropriately  combining  the  kernels. 
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Figure  4-7  is  the  ambiguity  diagram  for  the  kernel,  A^   and  B2, 
code  length  10,   (The  code  values  for  the  kernels  are  listed  at  the  end 
of  Appendix  F.)   The  ambiguity  diagram  is  the  term  used  to  describe  the 
three  dimensional  figure  generated  from  the  ambiguity  function.   The 
bandwidth*  of  the  ambiguity  function  in  Figure  4-7  is  approximately 
O.lir  .   For  the  system  postulated  in  Figure  2-2,  this  bandwidth  corresponds 
to  a  velocity  discrimination  of  3000  knots,  for  T  ■  1 u.sec,  fQ  =  10  GHz. 
For  fQ  =  1  GHz  and  T  =  0,1  usee,  velocity  discrimination  is  roughly 
300,000  knots.   Hence  complementary  codes  of  this  length  cannot  be  used 
for  situations  requiring  good  velocity  discrimination. 

Figure  4-8  is  the  ambiguity  diagram  for  the  kernel,   A2  and  B2, 
code  length  10,  and  with  5><j  in  increments  of  0.125"U"  . 

Figure  4-9  is  the  ambiguity  diagram  for  the  kernel,   A-i  and  Bj, 
code  length  10.   Note  that  the  range  sidelobes  extend  out  to  seven  units, 
whereas  the  range  sidelobes  in  Figure  4-8  extend  to  10  units.   However, 
the  amplitudes  of  the  sidelobes  for  the  kernel,  A^  and  B^  are  slightly 
larger  than  for  sequences  A2  and  B2. 

Figure  4-10  is  the  ambiguity  diagram  for  the  kernel,  A  and  B,  length 
26.   The  range  sidelobes  extend  to  14  units.   The  maximum  amplitude  of 
the  range  sidelobes  is  0.32,   The  ambiguity  diagrams  in  Figure  4-9  and 
4-10  are  in  the  same  increments  of   "(^  .   Note  that  the  bandwidth  of 
the  ambiguity  function  for  the  kernel,  A  and  B,  length  26,  is  less  than 
the  bandwidth  of  the  ambiguity  function  for  the  kernel,  A-,  and  Ei  length 
10.   This  effect  is  the  compression  phenomena  for  longer  codes  which  is 
discussed  in  Chapter  IV, 


*The  bandwidth  of  the  ambiguity  function  is  here  defined  as  doppler 
phase  (converted  to  frequency)  at  the  first  zero  crossing  of  the  ambiguitv 
function  for  T  =  0. 
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TABLE  1 

CONVERSION  FROM  DOPPLER  PHASE,   $j  ,  TO 
TARGET'S  RADIAL  COMPONENT  OF  VELOCITY,  vr. 


Note:   c  =  velocity  of  propagation  of  the  wave  in  the  medium 
(5.91  x  108  knots) 

f  =  carrier  frequency 

T  =  duration  of  a  code  symbol 

fd  =  doppler  frequency 

1\  =    TTfdT   =   2vrf0Ttr/C 


vr  (knots)       f   (GHz) 

.....  t   .,  .. 

T  (  fix  sec.) 

<£  d  (Trad) 

j 

60          1 

] 

. 

2.00xl0~4 

100 

3.38xl0"4 

600 

2.00xl0"3 

1000   j 

3.38xl0"3 

1600 

5.40xl0"3 

2000 

6.75x10-3 

3000   { 

1.00x10-2 

4000 

1.35x10";? 

5000 

1.69xl0~2 

6000 

t 

r 

i 

r 

2.00xl0-2 

60000 

\t 

, 

0.2 

60 

1 

0.1 

2.00xl0-5 

100 

3.38x10-5 

600 

2.00xl0"4 

1000 

3.38xl0"4 

1600 

5.40xl0~4 

2000 

6.75xl0~4 

3000 

l.OOxlO"3 

4000 

1.35xl0-3 

5000 

1.69xl0-3 

6000 

1 

r 

i 

r 

2.00x10-3 
2.00xl0-2 

60000 

f 

% 
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TABLE    1 
(cont'd) 


v   (knots) 

fo  (GHz) 



T  (  fj  sec.) 

^d  (Trad) 

60 

10 

1 

2.00xl0-3 

100 

3.38x10-3 

600 

2.00xl0-2 

1000 

3.38xl0-2 

1600 

5.40xl0"2 

2000 

6.75xl0~2 

3000 

0.1 

4000 

0.135 

5000 

0.169 

6000 
60000 

1 

0.2 

2.0 

60 

10 

0.1 

2.00xl0"A 

100 

3.38xl0"4 

600 

2.00xl0"3 

1000 

3.38x10-3 

1600 

5.40x10-3 

2000 

6.75x10-3 

3000 

1.00x10-2 

4000 

1.35x10-2 

5000 

1.69xl0-2 

6000 

2.00x10-2 

60000 

1 

f 

i 

' 

0.2 
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CHAPTER  V 

INVESTIGATION  OF  THE  DOPPLER  EFFECT 
ON  THE  TRANSMITTED  SIGNAL  (WIDEBAND  ANALYSIS) 


A.  BACKGROUND 

In  Chapters  II,  III,  and  IV  the  analysis  was  done  on  pulse-coded 
systems  for  the  narrow-band  case.   What  is  meant  by  narrow-band  is  that 
doppler  has  negligible  effect  on  the  envelope  of  the  transmitted  signal, 
but  shifts  the  carrier  frequency  by  the  doppler  frequency.   It  is  the 
purpose  of  this  chapter  to  discuss  the  doppler  effect  on  the  envelope 
of  the  waveform  and  to  derive  the  expression  for  the  time  waveform  of 
the  received  signal  (wideband  case).   As  stated  in  the  introduction, 
target  acceleration  is  assumed  to  be  negligible  over  the  interval  of 
measurement. 

B.  EFFECT  OF  DOPPLER  ON  THE  TRANSMITTED  SIGNAL'S  FREQUENCY  COMPONENTS 

Assume  that  s(t)  is  any  arbitrary  transmitted  signal.   Then  the 
Fourier  transform  (if  it  exists)  is  a  measure  of  the  frequency  content 
of    the  signal.   It  is  assumed  that  s(t)  is  sufficiently  well  behaved 
to       <  Fourier  transform.   If  the  shorthand  notation  is  utilized,  the 
Fourier  transform  of  s(t)  is 


J  [a(ol  =  £W  . 


It  is  well  known  that  the  doppler  effect  is  to  shift  the  frequency 
components  of  the  signal  such  that  the  new  frequency  f  is  given  by 

1  -V  Vr/C  ' 


r  =  (ii^M  4  5  5-1 


'•This  equation  is  a  statement  of  the  two-way  doppler  effect  for 
radar  and  sonar.   For  the  radar  case  Torres^-  has  indicated  that  if  the 
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where  v  =  radial  component  of  velocity  between  the  target  and  the 
transmitter 

c  =  speed  of  propagation  of  the  wave  in  the  medium 
Note:   The  convention  is  chosen  so  that  v  is  positive  for  closing 

targets  and  negative  for  opening  targets. 

Define  a  quantity,  the  doppler  factor,  D,  such  that 

f •  -  Df  t  5-2a 

Alternatively, 


0  =  i±Vg  5-2b 

l-Vc  • 


The  doppler  factor  may  be  given  in  terms  of  the  doppler  frequency  f  ,. 
Note  that  the  doppler  frequency  is  defined  as  the  change  in  the  carrier 
frequency  due  to  the  doppler  effect.   Hence 

fd  -  V  -  fo  .  5-3b 

Then 

fd  =  Dfo  "  fo  "  <D  "  Dfo  •  5"3c 


velocity  vectors  of  the  transmitter  and  the  target  are  noncollinear,  the 
angle  at  which  the  energy  arrives  at  the  receiver  (assuming  that  the  same 
antenna  is  used  for  transmitting  and  receiving)  differs  from  the  angle 
at  which  the  energy  was  transmitted  due  to  the  angular  aberration  phe- 
nomenon. 1" 

In  the  sonar  case  equation  5-1  applies  if  the  medium  is  at  rest 
with  respect  to  the  transmitter. 


64 


The  above  equation  may  be  rearranged  so  that 


D  =  1  +  |4  5-4a 


If  the  ratio  -£-  «  1,  then  the  doppler  factor  may  be  approximated  by 


D  =  1  +  M-  5-4b 

c   • 


hence  the  doppler  frequency  is  approximately 

f   =  £*£  f  5-5 

rd    c   ro  *  5  J 

Next  let  us  examine  the  received  signal's  Fourier  transform,  S'(f). 
(The  prime  notation  is  used  to  indicate  that  the  received  signal  differs 
from  the  transmitted  signal,)   Since  the  effect  of  doppler  is  to  shift 
all  the  frequency  components  of  the  original  signal  by  a  multiplicative 
factor,  D,  the  received  signal's  transform  is  given  by 

S'(f)  =  S(f/D)   .*  5-6 

C.   TIME  WAVEFORM  OF  THE  DOPPLER  SHIFTED  SIGNAL 

It  was  established  in  the  previous  section  that  the  received  signal's 
Fourier  transform  was  related  to  the  transmitted  signal's  transform  by 
a  change  in  argument.   Hence  the  inverse  transform  of  the  received  sig- 
nal's Fourier  transform  will  yield  the  time  waveform  of  the  received  signal. 


*We  arrived  at  the  amplitude  factor  being  unity  from  considerations 
of  the  special  theory  of  relativity.   In  Appendix  A  we  proved  that  the 
transmitted  and  received  energy  is  dependent  on  the  relative  velocity 
between  the  transmitter  and  the  target.   This  result  was  also  obtained 
by  A.  W,  Rihaczek,  12  Xo  obtain  the  energy  of  the  received  signal,  we 
integrate  the  squared  magnitude  of  equation  5-6  over  all  frequencies. 
Thus  the  energy  amplitude  factor  is  introduced  via  the  change  in  fre- 
quency scale. 
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Since  the  time  waveform  and  its  Fourier  transform  form  a  one  to  one 
correspondence,  the  time  waveform  ohtained  by  the  inverse  transform  of 
the  received  signal's  Fourier  transform  will  be  unique.   If  the  inverse 
Fourier  transform  is  taken,  then 


•w-  j'Ud)!-  ^a$i)J"M  .  5-7 


By  a  change  in  variable 


D 


-r-f   , 


then  equation  5-7  becomes 


»W  -  of  Ktt')  e'^V  . 


5-8 


Equation  5-8  has  the  same  form  as  the  transmitted  signal  with  the  excep- 
tion of  the  amplitude  factor  D  and  the  change  in  the  time  variable  to 
Dt.   Therefore 


,-i 


tfft>-  Ds(DO«  J  [S(&]  5-9 


Equation  5-9  is  a  well  known  transform  pair  and  is  referred  to  by 

13 
R.  M.  Bracewell    as  the  "Similarity  Theorem."   Since  the  received  sig- 
nal is  of  the  same  functional  form  as  the  transmitted  signal  except  for 
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a   change    in  argument   and    in   amplitude,    s'(t)   may   be  written  as 

s'(to    =    DKCft)  cc6  u)aD%  5-10a 


U 


No 


=  TXX 2  X^jllD^ -T)  -  AtDi -  (i+l)T]}  cos  u)0Dt ,  5-10b 


te   that    the   unit   step    function    jUftt-iY)    may  be  written   as   ^[oH-i^i]    and 


iLLoit-i^)l=.  ia-'^) . 


Define 


*D  =      »/D.  5"n 


Then 


<A>D  =     U^Ci+4^     ^       tJ&a+U^     ,  5-12 


■fi 


Hence  equation  5-10b  becomes 


s'(-t)  =  Dot  f  >t  { 1(1- LTD)  -ilft-U*0To]  ccsK+^  g    5-13 


i=i 


Equation  5-13  is  the  time  waveform  for  the  wideband  doppler  shifted 
signal.   Note  that  the  wideband  signal  reduces  to  the  narrow-band  signal 
if  the  doppler  factor  is  nearly  unity. 
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If  equation  5-13  is  examined  carefully,  we  see  that  doppler  produces 
three  effects. 

1)  The  carrier  frequency  is  shifted  by  the  doppler  frequency. 

2)  The  duration  of  a  code  symbol  is  altered  by  the  doppler  factor. 

3)  The  received  energy  differs  from  the  transmitted  energy  by  the 
doppler  factor. 

The  reason  for  the  carrier  frequency  being  shifted  by  doppler  is  obvious, 
and  Appendix  A  shows  the  relationship  between  the  received  and  transmitted 
energy.   However,  why  is  the  duration  of  a  code  symbol  altered  by  the 
doppler  effect? 

Let  us  determine  the  physical  significance  of  the  alteration  of  the 
duration  of  a  code  symbol. 
Let 

1  =  number  of  cycles  of  RF  in  the  duration  of  one  code  symbol 

of  the  transmitted  signal 
1    =  number  of  cycles  of  RF  in  the  duration  of  one  code  symbol 


Then 


of  the  received  signal. 


q    _  duration  of  one  code  symbol  for  the  transmitted  signal 

1    duration  of  one  cycle  of  the  transmitted  signal        5-14a 

=  T   _  f  T  5-lAb 

T7f;-f°T 

<l'    a   duration  of  one  code  symbol  for  the  received  signal    5-15a 
duration  of  one  cycle  of  the  received  signal 


=  f  'Tn  5-15b 

io  iD  • 


»V 


From  equation  5-11 


T 
T  =  — 
LD        D  , 
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and  from  equation  5-3b 


f  '  =  Df_ 
o      °  • 


Therefore 


^   ■  <Dfo> <£>  =  foT  ■  5"16 


Hence 


1  M   .  5-17 

Thus  we  see  that  the  number  of  cycles  of  RF  in  the  transmitted  and 
received  signals  remains  constant.   Kence  the  implication  of  the  doppler 
effect  is  that  when  the  transmitted  signal's  carrier  undergoes  a  shift 
due  to  doppler,  the  duration  of  a  code  symbol  changes  so  that  the  total 
number  of  cycles  of  RF  for  the  transmitted  and  received  signals  rerains 
constant.   This  is  compatible  with  the  special  theory  of  relativity  as 
shown  in  Appendix  D. 
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CHAPTER  VI 

WIDEBAND  AMBIGUITY  FUNCTION  FOR 
COMPLEMENTARY  CODES 


A.   BACKGROUND 

It  might  be  well  to  indicate  those  situations  in  which  the  doppler 
effect  on  the  envelope  of  the  transmitted  signal  cannot  be  ignored. 
Certainly  for  relative  velocities,  between  the  target  and  the  transmitter 
near  the  velocity  of  propagation  of  the  wave  in  the  medium,  this  effect 
cannot  be  ignored.   Another  situation  involves  the  duration  of  the  pulse 
train.   As  indicated  in  Chapter  V,  the  doppler  shifted  waveform  undergoes 
a  compression  or  expansion;  hence  if  the  total  duration  of  the  shifted 
waveform  is  compressed  or  expanded,  approaching  the  duration  of  one  code 
symbol,  energy  will  be  lost  in  the  correlation  process.   The  criterion 
for  the  use  of  the  narrow-band  analysis  is  given  in  Appendix  B  and  is 
duplicated  here.   Hence  if 


i-W  ,,   i 


«  IT  6-1 


■?  K       tf 


Or  alternatively, 


c    * 


NTW|vr\  ^<  -   ,     (See  equation  B-14)  6~;' 


Then  the  narrow-band  analysis  applies. 


*This  result  is  indicated  by  A.  W.  Rihaczek1  ,  and  E.  J.  Kelly  and 
R.  P.  Wishner.l* 
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Equation  6-2  is  a  statement  of  the  velocity-time-bandwidth-product 
uncertainty  relationship. 

B.  SYSTEM  BLOCK  DIAGRAM 

The  system  block  diagram  for  the  wideband  case  is  essentially  the 
same  as  the  narrow-band  case  with  the  exception  of  two  critical  areas. 
A  guard  interval  is  necessary  between  code  A  and  code  B  to  allow  for 
time  expansion  or  compression  of  the  two  codes  due  to  the  doppler  effect 
on  the  envelope  of  the  received  signal.   Also  the  delayed  replica  of  the 
transmitted  signal  must  be  clocked  at  the  rate  corresponding  to  the 
expected  clock  rate  of  the  received  signal  and  the  carrier  frequency  must 
be  shifted  by  the  doppler  frequency.   Suppose  that  the  replica  signal 
were  shifted  by  the  doppler  frequency,  but  not  clocked  at  the  expected 
rate;  then  degradation  of  the  output  of  the  matched  filter  would  be 
expected.   In  fact,  the  above  statement  implies  that  we  are  requiring  a 
narrow-band  system  to  detect  a  wideband  signal.   From  the  point  of  view 
of  equipment  complexity,  it  is  more  economical  just  to  shift  the  carrier 
frequency  of  the  transmitted  signal  to  obtain  the  replica  signal.   Appendix 
C  contains  the  analysis  of  the  replica  signal  just  being  shifted  in 
carrier  frequency  for  the  wideband  case.   It  is  expected  that  no  serious 
degradation  of  the  output  of  the  matched  filter,  and  hence  the  "envelope 
detector,"  will  be  evidenced,  so  long  as  the  time-bandwidth-product-velocity 
uncertainty  relationship  is  not  violated.   (Equation  6-2) 

C.  WIDEBAND  AMBIGUITY  FUNCTION  FOR  SINGLE  CODES 

In  the  wideband  case  the  replica  signal  is  obtained  by  clocking 
the  transmitted  signal  at  the  expected  rate  and  translating  this  signal 
by  the  expected  doppler  frequency.   For  convenience  the  wideband  matched 
filter  block  diagram  is  given  below. 
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s(t) 

i 

clocking  device  and 
frequency  shifter 


s'(t  -  tf1) 


^Si(Dit  -  rk) 

oh — — 


y^t  ^diAi') 


integrator 


■to  collator 


Wideband  matched  filter 
Figure  6-1 


The  replica  signal  is  given  by 


SM^  -  rk)  -  l(Dit;-rk)  cos  o)c(Dil-  rk)  , 


6-3a 


where  D . 

x 


expected  doppler  factor 
f. 


1  + 


-di 
fo 


*^di  =  expected  doppler  radian  frequency. 


The  received  signal  is  given  by 


s'(H*)=  D'xIdU-^)]  costK+u^M-W]  > 


where  D'  =  target's  doppler  factor 
fd' 


=  1  + 


fo 


'Ai'  =  target's  doppler  radian  frequency t 


6-3b 


6-4a 


6-4b 
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Then  the  output  from  the  matched  filter  is  given  by 

CO 

J  -oo 

=  D'  fltD'a-ta)]  cos  [( ^  utf  )(i-iA)}  o 

Now  let  us  shift  the  time  reference  by  letting 


tl'  t'ti    ' 


Thus  equation  6-5a  becomes 


Then  define 


6-5a 


<^lto*^lVV%-)l  d*i    .  6"3b 


1   Dj    d  • 


Since  t^  is  the  variable  of  integration,  it  may  be  replaced  by  t  with 
no  loss  in  generality.  Hence  the  output  of  the  matched  filter  may  be 
written  as 


oo 

-Oo 
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6-7a 


£-X   — Oft  -^ 


C  3 


+£•'  ^  ICo'tni^U  'q))  ooa[(jW^^-lv«^i)^J  &  0  6-7b 


The  second  term  in  equation  6-7b  will  be  rejected,  provided  that 
(l^j'  -  (J^ji)  <C<C  (2u50  +  U^')*   As  in  tne  narrow-band  case,  the  output 
of  the  matched  filter  is  an  ordinate  in  the  t    domain.   Then,  the  collator 
collects  these  ordinates  in  the  HT  domain.   The  output  of  the  collator 
is  given  by 


Let 


W0i  =  ^o+  %  6"9a 

^1=-  u^'-u^i   .  6-9b 

If  the  Euler  identity  for  the  cosine  function  is  used,  equation  6-8  becomes 

-co 
The  output  of  the  "envelope  detector"  will  be  given  by 


OO 
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To  the  extent  that  we  can  ascribe  an  envelope  to  the  output  of  the 
summer,  then 


?nv<^e{eiuVT'1|£iWt)i[o1U-r1^>^cl^)^ 


i^iA 


I  /OU  \   "til  L 

-OCi 


6-11 


In  Chapter  II,  it  was  proved  that  for  the  narrow-band  system  the  output 
of  the  "envelope  detector"  is  identical  to  the  ambiguity  function.  For 
the  wideband  case  the  ambiguity  function  will  be  defined  by 


-co 


6-12 


It  is  only  to  the  extent  that  the  approximation  made  in  equation  6-11 
holds,  that  it  is  possible  to  give  a  physical  interpretation  to  the 
wideband  ambiguity  function. 

Let  us  now  reference  the  received  signal  to  the  replica  signal  by 
letting 


H  ■  D*t. 


Thus  equation  6-12  becomes 


S^^di^a)  s 


-<x> 


6-13 


*If  the  integral  in  equation  6-11  is  a  "relatively"  slowly  varying 
function  of  Xx  »  as  compared  with  gjuWTi   (  then  this  approximation 
is  valid. 
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Then  define 


***** 


6-14 


r=  d^  v^a      , 


6-15 


AA  =  ^L     ^  j^'-  ^  \   g) 


/L0d-  i^jt  N 
Uo^  +  oIj.  / 


A'  — \  "-'o  > 


6-1.6 


=   tV^o  . 


6-16b 


Therefore  the  wideband  ambiguity  function  becomes 


,60  t   yVl 

A   -oo 


6-17 


Note  that  either  equation  6-17  or  equation  6-13  may  be  chosen  as  the 
definition  of  the  wideband  ambiguity  function. 

One  of  the  properties  of  the  narrow-band  ambiguity  function  is  that 
it  is  symmetrical  in  range  deviation,  t  ,  and  doppler,  iOj.   However, 
the  wideband  single  code  ambiguity  function  seems  to  exhibit  no  simple 
symmetry  property.   Thus,  we  will  obtain  equations  for  the  wideband 
single  code  ambiguity  function,  explicit  for  a  limited  case, to  indicate 
the  method  that  is  involved. 

The  wideband  single  code  ambiguity  function  will  be  calculated  for 
Uj^  >  0,  T>0.   These  equations  will  be  limited  to  the  case  of  the 
received  signal's  envelope  being  compressed  or  expanded  by  not  more  than 
the  duration  of  one  code  symbol  of  the  replica  signal.  Appendix  E 
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contains  the  analysis  of  the  limits  of  integration  and  the  code  values 
in  the  interval  of  integration  for  the  wideband t    single  code  amoiguity 
function.   We  have  divided  the  ambiguity  function  into  two  intervals  of 
range  deviation. 

1.   For  range  deviations,  *Z  t    in  the  interval 

we  have 


J]-t 


lri 


^Xl+J      ) 


U+fVr 


,i^ 


at 


t+(L-A)T0L 


U-lM 


1-1  Ci^)T' 


6-18 


where  N  =   length  of    the   code 


.  + 


^ 


X-  =  31  depending  on  the  ith  value  of  the  code  symbol 

D^  =  replica  signal's  doppler  factor  =  1  +  -^Li. 

f  o 

D'  =  received  signal's  doppler  factor  =  1  +  — ~ 

f  o 

D    =   PL  ■  f0  +  fd-' 
*     D±    fo  +  fdi 


T  =  transmitted  signal's  duration  of  a  code  symbol 


TDi 


T  f 


T'  = 


replica  signal's  duration  of  a  code  symbol  =  — 

Di 

received  signal's  duration  of  a  code  symbol  =  i 


V 


6-19 
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AT  =  magnitude  of  the  difference  in  the  duration  of  a  code  symbol 

T  .  -  T* I 
^Ui    l 


p  =  smallest  integer  >  Vt^  1 

j  =  J-l  (denotes  a  complex  quantity) 

0(  =  normalization  constant  =. 


"M 


NT  • 


-  6-19 
(cont.) 


Then  if  the  integration  is  performed,  equation  6-18  becomes 


i'Cr,^*  v^tf 


K-y 


X;X1+p  Uti^T^Ci-QTK-Tl 


i=l 


iiQ  ^rCi^rrUj-i^-Tl      5 ^ai+rirrfo^Hfl 


2. 


-uUnk-vq^vr'3 
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2.      For   range   deviations,   <C  ,    in   the    interval 


AT  AT         AT  "        AT 


m  =   1,2,",N   ; 


we   have 


5:'(?,uu 


DA0(' 


n-R-k 


IJUIOT 


^dl 


T+U-DTr 
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n-S-(ku) 
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^*sx 


i-L 


i-^1 


t+  irDj- 
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Cn-JklK' 


"*"  -^n-S-kX^.ji 


+1 


e1     d-b 


<>0t ' 


.N-* 


*r4(i-k>T5; 


+  2XA  e''dt 


i=M 


,i>U 


UHVr 


i=n-ft 


rUi-kH^ 
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where   the  definitions   in  equation   6-19  apply  and 


n   = 


AT 

1   AT 

q  =  smallest 

integer  > 

Vat 

1  =  q  modulo 

(n±) 

k  =  smallest 

integer  > 

^/tw 

If  the  integrals  in  equation  6-21  are  evaluated,  then  the  ambiguity 
function  becomes 


4'(r,u^  =  dacx2  I  2  ^iXi+k  Kuk^-q-i^-Ti 


sin  ^[(^k)Ti(q-irK.r3        .  ^[q.k^la-irpil 


2x 


eJ 


n-JWku) 


1=1. 
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Sin  ^  A    2,       e      ^ 


6-22 


D.   WIDEBAND  AMBIGUITY  FUNCTION  FOR  COMPLEMENTARY  CODES 

As  stated  in  Chapter  IV-C,  to  obtain  the  ambiguity  function  for  a 

complementary  code,  given  the  ambiguity  function  for  single  codes,  one 

merely  replaces  the  signed  product  for  the  single  codes  by  the  signed 

products  for  the  complementary  codes. 

1.   Consider  range  deviation,  *T  ,  in  the  interval 
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Then  if  the  signed  products  for  the  single  codes  in  equation  6-20  are 
replaced  by  the  signed  products  for  the  complementary  codes,  the  ambi- 
guity function  for  the  complementary  codes  will  be  given  by 


5>»»J  -  DAoC 
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[XA|XAw-v  X6iX6i+i> 
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where      q( 


c      >| 


JL 
NT  • 


2,   Now  consider  range  deviation,  T  ,  in  the  interval 


A  l  AT    AT     AT 


in  if  the  signed  products  in  equation  6-22  are  replaced  by  the  signed 
products  for  the  complementary  codes,  the  ambiguity  function  for  the 
".omplementary  codes  will  be  given  by 
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CHAPTER  VII 
CONCLUSIONS 

A.   SUMMARY 

The  research  presented  in  this  thesis  had  two  primary  objectives. 
The  first  objective  was  to  formulate  the  ambiguity  function  of  the  comple- 
mentary codes  for  the  narrow-band  case.   The  second  objective  was  to 
formulate  the  ambiguity  function  of  the  complementary  codes  for  the 
wideband  case.   We  shall  summarize  the  results  of  these  investigations, 
but  not  necessarily  in  the  order  presented  in  the  thesis. 

1.   NARROW-BAND  ANALYSIS 

In  the  narrow-band  case,  Siebert's  ambiguity  function  was  extended 
to  include  complementary  codes.   The  ambiguity  function  of  complementary 
codes  was  defined  so  that  its  autocorrelation  function  contained  no  range 
sidelobes.   From  the  definition  we  determined  that  the  two  codes  should 
be  combined  at  the  outputs  of  the  matched  filters  to  effect  the  necessary 
cancellation  of  range  sidelobes.   During  the  course  of  research  into  the 
ambiguity  function  of  complementary  codes,  we  discovered  that,  once  the 
ambiguity  function  for  a  single  codes  has  been  derived,  then  to  derive 
the  ambiguity  function  for  a  complementary  code  is  a  simple  task  —  replace 
the  signed  products  of  the  single  code  with  the  sum  of  the  signed  products 
of  sequences  A  and  B,  i.e., 

replace  (X^j)  by  (XA^XAj  +  XB^XBj )  . 

In  the  derivation  of  the  ambiguity  function  for  single  codes,  it 
was  proved  that  the  ambiguity  function  was  identical  to  the  output  of 
the  "envelope  detector."   (The  "envelope  detector"  refers  to  the  device 
in  the  system  that  takes  the  envelope  of  the  signal,  but  in  the  range 
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or  T  domain.   This  use  of  the  term  differs  from  conventional  terminology 
in  that,  usually,,  envelope  detection  is  done  in  the  time  domain.)   Although 
this  result  is  not  unique,  it  helps  to  clarify  some  of  the  misinterpre- 
tations of  the  ambiguity  function.   It  is  a  common  misconcepticr  to  state 
that  the  output  of  the  matched  filter  is  identical  to  the  ambipuity 
function. 

One  result  obtained  in  the  analysis  of  the  ambiguity  function  of 
complementary  codes  is  that  in  the  doppler  domain,  the  zero  crossings 
of  the  ambiguity  function  move  closer  to  the  origin  as  the  length  of  the 
code  increases,  if  the  duration  of  a  code  symbol  is  held  constant.   The 
implication  of  this  result  is  that  longer  codes  permit  better  velocity 
d^^rimination.   However,  this  process  of  increasing  the  code  length  to 
obtain  better  velocity  discrimination  cannot  be  continued  indefinitely 
without  incurring  a  penalty.   The  penalty  is  that  if  the  velocity-time- 
bandwidth-product  uncertainty  relationship  is  violated,  then  distortions 
of  the  envelope  of  the  received  signal  must  be  taken  into  account.   This 
amounts  to  having  to  clock  the  replica  signal  at  the  expected  rate  and 
hence  increasing  the  complexity  of  the  system. 

2.   WIDEBAND  ANALYSIS 

In  the  wideband  analysis,  it  can  be  shown  that  the  doppler  effect 
is  three-fold. 

a.  The  carrier  frequency  is  shifted  by  the  doppler  frequency. 

b.  The  envelope  of  the  transmitted  signal  is  compressed  or 
expanded  according  to  whether  the  target  is  closing  or 
opening  with  respect  to  the  receiver. 

c.  The  received  energy  differs  from  the  transmitted  energy 

by  the  same  factor  as  the  received  and  transmitted  carrier 
frequencies. 
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Although  all  three  effects  can  be  proved  from  the  special  theory  of 
relativity,  only  the  last  two  were  proved  in  this  manner,   The  distortion 
of  the  envelope  of  the  received  signal  was  also  proved  as  a  logical  con- 
sequence of  the  application  of  Fourier  analysis,  and  the  doppler  effect 
on  the  carrier  frequency.   A,  W.  Rihaczek12  proved,  by  considering  a 
photon,  that  the  received  and  transmitted  energies  differ.   Our  approach 
to  the  problem  was  to  make  a  Lorentz  transformation  of  the  field  tensor. 
It  was  proved  that  for  binary  phase-reversal  modulated  signals,  if 

NTW|vr|  «   |  , 

then  the  narrow-band  analysis  applies. 

The  itfideband  ambiguity  function  was  derived  here;  however  this 
result  had  been  previously  accomplished  by  A.  W.  Rihazek^,  and  Kelly 
and  Wishner   .   We  merely  clarified  certain  aspects  of  the  analysis. 
Then  the  equations  for  the  single  code  and  complementary  code  wideband 
ambiguity  functions  were  derived.   These  equations  are  limited  to  the 
case  of  signal's  envelope  being  altered  by  not  more  than  the  duration 
of  one  code  symbol.   They  are,  indeed,  formidable,  and  a  computer  is 
required  to  obtain  the  ambiguity  diagrams.   Due  to  time  limitations,  the 
wideband  ambiguity  diagrams  were  not  completed. 
B.   AREAS  OF  POTENTIAL  STUDY 

Many  areas  of  future  research  folloxj  directly  from  this  thesis. 
A  logical  extension  of  this  thesis  would  be  to  construct  the  wideband 
ambiguity  diagrams  for  complementary  codes  and  single  codes  (Barker  codes, 
truncated  pseudo-noise  codes,  etc.)   Also  the  region  of  validity  for  the 
narrow-band  analysis  should  be  more  precisely  defined.   We  found  that, 
in  the  preliminary  analysis  of  the  computer  runs  for  the  wideband  ambi- 
guity function  of  complementary  codes,  the  wideband  ambiguity  function 
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for  complementary'  codes  began  to  deviate  substantially  from  the  narrow- 
band ambiguity  function  at  approximately 

NTwIvl  =  0.1  c/2 
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APPENDIX  A 

ENERGY  OF  THE  RECEIVED  SIGNAL 

DERIVED  FROM  CONSIDERATIONS 

OF  THE  SPECIAL  THEORY  OF  RELATIVITY 


Consider  the  reference  frame,  S,  with  the  transmitter  located  at 
the  origin  of  the  coordinate  system.   Assume  that  the  target  is  a  perfect 
reflector,  which  is  moving  at  a  velocity,  vr,  relative  to  the  transmitter. 
The  target  is  situated  at  the  origin  of  the  reference  system,  S'.   A 
suitable  set  of  axes  is  chosen  so  that  the  x  direction  of  both  of  the 
systems,  S  and  S1,  is  collinear  with  the  relative  velocity  vector,  v  . 
Figure  A-l  reflects  a  suitable  choice  of  axes.   Note  that  in  Figure  A-l, 
the  convention  for  positive  velocity  is  away  from  the  transmitter  and  is 
a  standard  convention  for  most  physics  textbooks.   Although  this  con- 
vention is  contrary  to  the  one  chosen  in  the  main  text  of  this  thesis, 
no  confusion  should  arise  here. 

The  transmitter  emits  a  plane  wave  of  finite  duration,  NT.   Let  field 
quantities  as  seen  in  the  transmitter  reference  frame  be  Ey  and  H?. 

Next,  we  require  the  field  quantities  as  seen  by  the  target  in 
reference  system,  S'„   These  quantities  can  be  obtained  by  application 
of  the  Lorentz  transformations  of  the  field-strength  tensor-"--*.   They  are 
given  by 

V  =  Ex  "  °  V  ■  Hx  -  ° 

V  =  *<Ey  -  £HZ)  Hy'  -  tf(Hy  +  £EZ)  =  0     A-l 

V  =  ^(Ez  +  £Hy>  =  °         V  "  3(Hz  "  $V  » 
where 


&=  Vr/c 


■fc   .-   **-& 
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transmitter 


S' 

y' 

target 


H, 


vr 


Reference  frames  for  the  target  and  the  transmitter 

Figure  A-1 


S 

y 


S* 

y' 


The  plane  wave  after  reflection 
Figure  A-2 


Hz  t  t  t 


Eyl  M 


Q  t 


yT 


The  plane  wave  after  reflection 
as  viewed  by  the  transmitter 

Figure  A-3 
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Then  we  want  the  reflected  wave's  field  components  as  viewed  by 
the  target.   These  field  components  can  be  easily  obtained  by  applying 
the  boundary  conditions  for  a  perfect  reflector.   Therefore,  the  field 
quantities  as  seen  by  the  target  are  given  by 

E  "  -  0  II  ''  =  0 

x  x 

E  "  =  -F  '  tL."  =0  A-4 

y      y  y 

E  "  =  0  H  '  '  =  +11  '  . 

z  z       z 

Finally,  we  require  the  field  quantities  as  viewed  by  the  transmitter. 
The  reflected  wave's  components  can  be  easily  obtained  by  applying  the 
transformation  rules  in  equation  1,  and  noting  that  the  velocity  component 
is  now  in  the  negative  x  direction;  hence  we  must  replace  vr  by  -vr. 
Let  us  denote  the  reflected  field  quantities  by  a  triple  prime  notation. 
Thus 

L  '"  =  0  1^"'  =  0 

L/"  -  *(Ey"  +  eV»)       y '  -  0  A-5 

Ez'"  =  0  Hz'"  =  ^(Ha"  +  My"). 

If  equation  A-4  is  substituted  into  equation  A-5^ 

Ey'"    =     ^("Ey'    +    ^llz») 

liz'"  =  S(HZ'  -  f>Ey'),  A-6 

Then   substituting   equation  A-l    into   equation  A-6 

E    '"    =       *{[-tf(Ey   -     f>HzJ[     +        ^(llz   -    $ly)\ 


y 

=  _  v2n    +    A2^ 


=  -  ^2(1   +    $2)E     +  2  ^2^>EZ  A-7 


Hz'"    =        *{*(HZ    -     £Ey)-       ^(Ey   -    &H,)} 
=      o2U   +     ^>2)HZ   -   2tf2£,Ey 
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A-8 


The  solutions  for  Maxwell's  equations  are  such  thac 


E  -  2l  -  H   • 


since 


Z  =   -5s     ■  1  ,  if  the  field  quantities  are  in  Gaussian  units. 

O     -\J  t0 


Therefore,  equation  A-7  becomes 


.'"  -  L-^2(l  +  £2)  +  2  *2£>]  E, 

=  -x2u  -  e>)2E   , 


A-9 


and 


Hz"»    =     *2(1  +    ^2)HZ   -   2  S2£K. 
=     tf2U   "   £)2HZ    . 


A-10 


The  power  density  can  be  obtained  by  Poynting's  Theorem,  and  is 
given  by 


p«-k&{i^*H'i  , 


A-ll 


Thus  the  power  density  of  the  reflected  wave  as  seen  by  the  transmitter  is 


A- 12 


e  E   is  the  amplitude  of  E  . 

o  '  y 
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The  power  density  of  the  transmitted  wave  is 


p  I  =  -M  e\z  a"13 


8tt 

The  energy  density  of  the  wave  can  be  determined  by  integrating 
the  power  density  over  the  duration  of  the  signal.   Note  that  the  power 
density  is  uniform  over  the  planar  surface.   Thus  the  energy  density  of 
the  transmitted  wave  is 


go  2 

$  =    [   \9\dh  =  fJzQ\  NT.  A-L4 


And  the  energy  density  of  the  reflected  wave  is  given  by 


OO  1 


where  N  ■  length  of  the  code 

T,,f  =  duration  of  a  code  symbol  for  the  received  signal 

T  =  duration  of  a  code  symbol  for  the  transmitted  signal. 

During  the  analysis  presented  in  Chapter  V  and  Appendix  D,  it  is  shown 
that  the  received  signal's  duration  of  a  code  symbol  is  altered  by  the 
doppler  effect  so  that 

T 


T '  '  '  = 


D   ' 


where 


-  +  Yr/C  . 
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Thus   equation  A-15   becomes 


£'"--t^X^)zil) 


=  -^  la0\zDiiT    .  a-16 


Then  the  ratio  of  the  received  signal's  energy  density  to  the  transmitted 
signal's  energy  density  is  given  by 


_    =  JL    z 

%Tr\Eo\  U  I 


=  D  .  A-17 
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APPENDIX  13 

CRITERION  FOR  THE  USE  OF  THE 
NARROW-BAND  APPROXIMATION 
(VELOCITY-TIME- BANDWIDTH-PRODUCT  UNCERTAINTY  RELATIONSHIP) 


In  this  appendix  we  shall  derive,  from  rather  elementary  consider- 
ations, the  region  of  validity  for  the  narrow-hand  approximation.   This 
approximation  is  made  in  Chapters  II  and  III  and  assumes  that  doppler 
has  negligible  effect  on  the  envelope  of  the  signal.   Chapter  V  contains 
the  derivation  of  the  time  waveform  of  the  doppler  shifted  signal  for 
the  wideband  case.   (For  purposes  of  clarification,  the  terms,  doppler 
shifted  signal  and  the  received  signal,  are  used  synonymously.)   We  con- 
cluded that  the  effect  of  doppler  was  threefold: 

1.  The  duration,  T^,  of  a  code  symbol  for  the  doppler  shifted 
waveform  was  altered  so  that 

TD  =  T/D.  B-l 

2.  The  carrier  frequency  was  shifted  by  the  doppler  frequency. 

3.  The  received  signal's  energy  density  differed  from  the  trans- 
mitted signal's  energy  density  by  the  doppler  factor. 

If  the  doppler  shifted  signal's  duration  of  a  code  symbol  is  altered  by 
doppler  so  that  the  total  duration  of  the  signal  is  expanded  or  compressed 
by  the  duration  of  a  code  symbol  of  the  transmitted  signal,  we  expect 
degradation  of  the  output  of  the  matched  filter.   In  particular,  if  the 
code  is  completely  random,  then  for  zero  range  deviation  the  output  of 
the  matched  filter  will  be  zero.   Hence  we  shall  use  as  our  criterion  of 
the  narrow-band  approximation,  that  the  total  duration  of  the  received 
signal  must  be  compressed  or  expanded  by  much  less  than  the  duration  of 
a  code  symbol  of  the  transmitted  signal,  T.   Let  us  examine  two  cases. 
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1,   Case  I:   The  total  duration  of  the  received  is  compressed  by 
much  less  than  the  duration  of  a  code  symbol  of  the  transmitted  signal. 
(This  statement  implies  that  positive  doppler  is  present.) 

Figure  B-l  indicates  the  envelopes  of  the  transmitted  signal,  I(t), 
and  the  received  signal,  I(Dt),  for  some  arbitrary  code.   lence  for  the 
condition  that  the  total  duration  of  the  received  signal  is  compressed 
by  much  less  than  the  duration  of  a  code  symbol  of  the  transmitted  signal, 
the  following  inequality  applies  : 

NTD  »  (N  -  1)T  .  E-2 

If  the  above  inequality  is  rearranged, 

T/TD  ^  (N/(N  -  1)  .  E-3 


From   equation   B-l 


Therefore 


And 


T/TD   =  D     «  K/(N  -   1)    , 


1   +  fd/fo     ^    N'/(N  ~   D    • 


fd/fQ       «    1/(N  -   1)    .  B-4 


Alternatively, 


1  +  v    /c 

- —     «    K/(N  -   1)    . 

1   -  vr/c 


If   N    »    1,    then 


vr/c      «  1/2K    .  E-5 
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Kote  that  the  time-bandwidth  product  is  equal  to  N»* 
Hence 

NTW  v   «.  c/2  .  B-6 

2.   Case  II:   The  total  duration  of  the  received  signal  is  expanded 
by  much  less  than  the  duration  of  a  code  symbol  of  the  transmitted  signal. 
(This  statement  implies  that  negative  doppler  is  present.) 

Figure  B-2  indicates  the  envelopes  of  the  transmitted  signal  and 
the  received  signal  for  some  arbitrary  code.   Hence  for  the  condition 
that  the  total  duration  of  the  received  signal  is  expanded  by  much  less 
than  the  duration  of  a  code  symbol  of  the  transmitted  signal,  the  following 
inequality  applies: 

NT  »  (K  -  1)TD  .  B-7 

The  above  inequality  may  be  rearranged  so  that 

T/TD  »  (N  -  1)/N  .  B-8 

Note  that  the  doppler  frequency  is  negative  in  this  case.   Fence  the 
doppler  factor  may  be  written  as 

\fd\ 

D  =  1 .  B-9 

fo 


*V'e  have  defined  bandwidth  as  the  video  bandwidth  to  the  first  null 
of  the  rectangular  pulse  spectrum,  i.e.,  W  =  1/T.   This  is  approximately 
equal  to  the  "effective"  RF  bandwidth.     Thus  NTW  =  N.   The  time-band- 
width product,  as  used  by  Torres,  was  defined  in  such  a  fashion  (using 
standard  deviation  in  both  time  and  frequency  domains)  as  to  be  a  complete 
and  precise  description  for  a  Gaussian  pulse.   If  rectangular  or  other 
shaped  pulses  are  used,  a  numerical  factor  of  the  order  of  magnitude  of 
unity  may  be  introduced  into  his  precise  inequality.   Hence  there  exists 
an  apparently  discrepancy  of  a  factor  of  2  between  Torres'  inequality 
and  our  inequality  B-6. 
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If  equation  B-9  is  substituted  into  the  inequality  B-P,  then 

fo        N 
f    ^   M  «  B-10 


Alternatively,  the  doppler  factor  may  be  given  by 

jvrj 

D  =    ■     "   C 

1+ .  JW  .  B-ll 

c 

hence,  substituting  equation  B-ll  into  the  inequality  B-8,  then 


\vrl  ^   1 

-c —  *  2N  -  1  .  E-12 


If    the   length   of    the   code   is  much   greater   than  unity 


w£L«1 


2N      .  B-13 

Then,  we  note  that  the  time-bandwidth  product  is  equal  to  the  length 
of  the  code,  N.   Hence 

NTW  \vr\  «  |  .  B-14 

Thus  if  the  time-bandwidth  product  multiplied  by  the  radial  component 
of  relative  velocity  is  much  less  than  one-half  the  velocity  of  the  wave 
in  the  medium,  then  the  narrow-band  assumption  applies.   Inequalitv  B-14 
is  a  statement  of  the  velocity-time-bandwidth-product  uncertainty  relation- 
ship. 
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APPENDIX  C 

EFFECT  OF  DOPPLER  MISMATCH,  DUE  TO  THE  REPLICA 
SIGNAL  ONLY  BEING  TRANSLATED  IN  FREQUENCY,  ON 
THE  OUTPUT  OF  THE  "ENVELOPE  DETECTOR" 


1.  BACKGROUND 

For  the  wideband  case,  the  system  block  diagram  contained  range 
banks,  such  that  in  one  range  bank,  the  received  signal  is  correlated 
with  an  exact  replica  of  itself.   Note  that  the  replica  signal  has  the 
same  carrier  frequency  and  clock  period  as  the  received  signal.   Hence, 
in  generating  this  signal  from  the  transmitted  signal,  the  system  must 
contain  a  clocking  device,  such  as  a  shift  register,  as  well  as  a  fre- 
quency translator.   Naturally  the  clocking  device  adds  to  the  complexity 
and  hence  the  cost  of  the  system.   Suppose  that  in  trying  to  reduce  the 
complexity  of  this  system,  we  merely  translate  the  frequency  of  the 
transmitted  signal  by  the  doppler  frequency.   How  is  the  output  of  the 
"envelope  detector"  affected?   If  the  replica  signal  is  obtained  by  a 
mere  frequency  translation,  we  are,  in  effect,  requiring  a  narrow-band 
system  to  detect  a  wideband  signal.   Thus  we  expect  no  serious  degradation 
of  the  output  of  the  envelope  detector,  so  long  as  the  time-bandwidth- 
product  uncertainty  relationship  (equation  6-2)  is  not  violated. 

2.  OUTPUT  OF  THE  ENVELOPE  DETECTOR  (SINGLE  CODE  SYSTEM) 

We  will  investigate  the  single  code  system  prior  to  the  complementary 
coded  system,  since  once  the  equations  for  the  single  code  system  are 
obtained,  it  is  a  simple  task  to  obtain  the  equations  for  the  complementary 
coded  system.   Figures  2-2a  and  2-2b  will  be  used  as  the  system  block 
diagram,  since  our  receiver  is  a  narrow-band  system.   The  received  signal, 
s'(t),  will  be  given  by  equations  5-10a  and  5-10b.   The  replica  signal, 
sM(t),  is  simply  the  transmitted  signal  translated  by  the  doppler  frequency} 
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hence  it  is  identical  to  the  narrow-hand  received  signal.   Then  s>j(t) 
is  given  by  equation  2-4.   Thus 


s'Cfc)  =   D'l(O'-fc)  Cos(i0o4-uud')\,  ? 


C-la 


s^&w  iauosdoe+^t  , 


C-lb 


Then  the  output  of  the  mismatched  filter  in  Figure  2-2b  is  given  by 


V^^i^cO 


oo 


-  _D_  tOS  (00T  ^  iCD'-feiJtt-T)^ 
2k       -0(J 


«to-«cy 


>  ,00 


+  5- \  KO'^ICV^  Cc^C^^^-C^^t]^  .        C-2 


-oo 


The  second  term  in  equation  C-2  will  be  rejected  by  the  integrator. 
Then  the  mismatched  output  of  the  "envelope  detector"  will  be  given  by 


XM(r,tf>at^') 


"M 


C-3 


%-«aa' 


From  equation  6-17  we  see  that  the  mismatched  output  from  the  "envelope 
detector"  is  related  to  the  wide-band  ambiguity  function  by 


X^fe,^*  $'(r>«o 


C-4 


D>1 
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Recall  from  equation  6-14  that 


DA  =  D'/Di# 


Then  the  differential  doppler  factor,  in  this  case,  is 


C-5 


Next,  let  us  obtain  the  mismatched  output  of  the  "envelope  detector" 
for  the  case  of  the  replica  signal's  doppler  frequency  matched  to  the 
received  signal's  doppler  frequency.   This  can  be  obtained  by  setting 
D.*  =  1,  and    <*)A  =  ^»  ^n  equations  6-20  and  6-22. 

1.   For  range  deviations,  T  ,  in  the  interval 


(m- 


-^^It^-^Tt  ; 


m  =  1,2,' ••,K> 


we  have 


M- 


A 


C-6 


where  N  =  length  of  the  code 

X^  =  tl   depending  on  the  ith  value  of  the  code  symbol 
D'  =  received  signal's  doppler  factor 


-  1  + 


D.  =  D' 


Lll 


r 


C-7 
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T  =  transmitted  signal's  duration  of  a  code  symbol 
Tn'  =  received  signal's  duration  of  a  code  symbol 

=  T/D' 
AT  =  magnitude  of  the  difference  in  the  duration  of  a  code 
symbol  -   |T  -  TD'\ 
p  =  smallest  integer  >\  yyy    , 


r 


C-7 
(cont. ) 


2.      For  range  deviations,    t  ,    in  the   interval 


We  have 


li--^-vl-v(m-l)Z  ^  H.  L  m^.      . 
AT  At     at         at      * 


m  =  1,2,  ••",::. 


n-S-k 


X^CJi  =   00(2|  Tl  X:X;.^[(Uk)^o-(i-l)T-r] 

z 


Ul 


tvA-(kn) 

+2  XiXUk+1  Er+vr-(utfT-p'] 


N-z 


+2  Xi.^4..  [r*U-k)r- (i+ifV] 


1*«-A 


tf-Z 


Ut\4 


C-8 
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where  the  definitions  in  equation  C-7  apply  and 


n  = 
nl   = 

q  = 
k  = 


TD7  AT 


T/A.T 

smallest   integer  >  \t/^r\ 

smallest    integer  >  \*r/T\    , 


J 


C-9 


3.   OUTPUT  OF  THE  ENVELOPE  DETECTOR  (COMPLEMENTARY  CODED  SYSTEM) 

As  indicated  in  Chapter  4-C,  it  is  easy  to  ohtain  the  equations  for 
the  complementary  coded  system  once  the  equations  for  the  single  code 
system  are  obtained  by  simply  replacing  the  signed  products  of  the  single 
codes  by  the  signed  product  for  the  complementary  codes.   Hence  if 

X^Xj  is  replaced  by  XA-jXAj  +  XE^XBj 

in  equations  C-6  and  C-8;  then  the  resulting  equations  will  be  the  mis- 
matched output  of  the  envelope  detector  for  the  complementary  coded 
system. 


105 


APPENDIX  D 

DURATION  OF  A  CODE  SYMBOL  DERIVED  FROM  CONSIDERATIONS 
OF  THE  SPECIAL  THEORY  OF  RELATIVITY 


In  Chapter  V  it  was  proved  that  the  transmitted  signal's  duration 
of  a  code  symbol,  T,  differed  from  the  received  signal's  duration  of  a 
code  symbol,  T',  by  the  inverse  of  the  doppler  factor,  D.   This  result 
was  obtained  from  considerations  of  the  doppler  effect  on  the  carrier 
frequency  and  also  by  the  use  of  the  Fourier  transform. 

In  this  Appendix  we  will  prove  the  same  result,  but  from  consider- 
ations of  the  special  theory  of  relativity. 

Consider  the  reference  frame,  S,  with  the  transmitter  located  at 
the  origin  of  the  coordinate  system.   Next,  consider  a  moving  frame,  S', 
with  the  target  located  at  the  origin.   A  suitable  set  of  axes  is  chosen 
so  that  the  x  direction  of  both  S  and  S'  is  collinear  with  the  velocity 
vector,  vr.   At  some  initial  time,  i.e.,  t  =  0,  we  shall  assume  that  the 
two  frames  are  coincident.   At  t  =  0,  the  transmitter  emits  a  plane  wave 
pulse  of  duration,  T.   The  pulse  travels  at  the  speed  of  light.   Then  we 
can  view  the  pulse  as  extending  in  space  so  that  its  trailing  edge  has 
an  x  coordinate  of  -cT.   For  t  >  0,  the  moving  frame  travels  in  the  +x 
direction  with  a  velocity  v  ,   Note  that  the  convention  for  positive 
velocity  is  away  from  the  transmitter  and  is  a  standard  convention  for 
most  physics  textbooks.   Although  this  convention  is  contrary  to  the  one 
chosen  in  the  main  text  of  this  thesis,  no  confusion  should  arise  here. 

At  time  t  =  T,  the  trailing  edge  of  the  pulse  reaches  the  origin 
of  S,  and  the  leading  edge  of  the  pulse  is  now  at  x  =  cT.   Now  let  us 
consider  the  spacial  distribution  of  the  pulse  in  the  moving  frame.   In 
the  moving  frame  the  x  coordinate,  x' ,  and  the  time  in  this  frame,  t', 
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can  be  derived  from  their  counterparts  in  the  reference  frame  from  the 
Lorentz  transformations.  "  Thus 


.»  _ 


=  *(x  -  vrt) 

t'  =   *<t  -It  x), 

c^ 


D-l 
D-2 


where 


*=H?H 


-1/2 


Let  us  denote,  as  in  Figure  D-l,  the  trailing  edge  of  the  pulse  as  event  1, 
For  event  1,  the  x  coordinate  of  the  trailing  edge  of  the  pulse  in  the 
moving  frame  is  given  by 


t(x   -  vrt) 


*vrT  . 


D-3 


For  event  1,  the  time  elapsed  in  the  moving  frame  is  given  by 


cl  = 


=  tf(t-^x) 


=  *T 


D-4 


Since  the  pulse  is  traveling  at  the  speed  of  light,  the  total  time  for 
the  trailing  edge  of  the  pulse  to  pass  the  observer  at  the  origin  of  the 
moving  frame  is  given  by 


1   C  G      * 
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(1) 


x=0 
t=T 


x=cT 
t=T 


1) 


xi=-  *vrT  /x'=0 


y' 


Reference  frames  for  the 
transmitter  and  the  target 

Figure  D-1 


Reference  frames  for  the 
target  and  the  receiver 

Figure  D-2 
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1 


1  -  vr/c 


T.  D-5 


Equation  D-5  is  a  statement  of  the  expansion  or  compression  of  the 
duration  of  a  pulse  for  the  case  of  one-way  doppler. 

Now  we  would  like  to  consider  the  duration  of  the  pulse  which  is 
reflected  from  the  target,  as  seen  by  an  observer  (the  receiver)  at  the 
origin  of  the  reference  frame,  S.   To  avoid  confusion  let  us  denote  this 
reference  frame  as  S'';  however,  bear  in  mind  that  we  are  assuming  that 
the  transmitter  and  receiver  are  at  the  same  physical  location.   We  can 
view  the  reflected  wave  as  being  retransmitted  from  the  target. 

Consider  S' ,  the  target's  frame,  and  S'',  the  receiver's  frame,  to 
be  coincident  at  t'  =  0,   From  S',  the  receiver  seems  to  be  moving  at  a 
velocity  vr  in  the  negative  x  direction.   Then  at  t'  =0,  the  leading 
edge  of  the  pulse  is  at  the  origin  of  both  frames.   At  t'  =  T',  the 
trailing  edge  of  the  pulse  is  at  the  origin  of  the  frame,  S'.   Let  us 
denote  this  as  event  2.   To  obtain  the  x  coordinate  of  the  trailing  edge 
of  the  pulse  as  seen  by  the  observer  on  S'',  we  utilize  the  Lorentz  trans- 
formations as  given  by  equations  D-l  and  D-2,  but  remembering  to  replace 
v„  by  -vr.   Thus  for  event  2 


=   *dT\ 


x.  =o 


D-6 
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For  event  2   the  time  elapsed  in  S'1  is 


\£   =  *(V-  |x') 


V-T 


-  St 


D-7 


Then  the  total  time  required  for  the  trailing  edge  of  the  pulse  to  reach 
the  observer  on  S' '  is  given  by 


■II  4   "  .  >< 


H^  -  *r'+ 


«VfT' 


*(i*£yr' 


^ 


l+^r/c   -y' 


1  -  Vc 


D-8 


Then   if   equation  D-5    is   substituted    into   equation  D-8 


T 


in 


i  +  \(c 


T/D    , 


D-9 


where 


0  - 


_      1-Tf/c 


14Vr/c 
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The  factor  which  compresses  or  expands  the  duration  of  the  pulse  is 
identical  to  the  one  ohtained  by  the  Fourier  transform  analysis  done  in 
Chapter  V,  with  the  exception  of  the  difference  in  sign  of  the  velocity 
due  to  the  change  in  sign  convention  adapted  in  this  Appendix.   In  our 
initial  assumptions,  we  assumed  that  the  duration  of  the  transmitted 
pulse  was  T.   It  is  easy  to  see  that  if  a  series  of  N  pulses,  each  of 
duration  T,  were  transmitted,  then  the  compression  or  expansion  factor 
would  be  identical  to  the  one  in  equation  D-9.   Hence  we  have  proved  that 
the  doppler  effect  alters  the  duration  of  a  code  symbol  by  the  factor  D. 
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APPENDIX  E 

DETERMINATION  OF  THE  LIMITS  OF  INTEGRATION 
FOR  THE  WIDEBAND  SINGLE  CODE  AMBIGUITY  FUNCTION 


In  this  Appendix  we  would  like  to  determine  the  limits  of  integration 
and  their  corresponding  code  values  in  the  interval  of  integration  for 
the  wideband  single  code  ambiguity  function.   In  the  formulation  of  the 
equations  for  the  wideband  single  code  ambiguity  function,  there  are 
four  quadrants  to  consider.   We  will  consider  only  the  first  quadrant, 
i.e.,    t£A>0,   f  >0,  and  will  limit  the  validity  of  the  equations  to 
the  case  of  the  received  signal's  total  pulse  duration  being  expanded 
or  compressed  by  not  more  than  the  duration  of  one  code  symbol  of  the 
replica  signal. 

By  equation  6-17,  the  wideband  single  code  ambiguity  function  is 
defined  as 


it'foaj-  \  ^fi(oAOitt  W3^|  . 


z 


E-l 


Hence  to  determine  the  limits  of  integration  and  the  code  values  in  the 
interval  of  integration,  we  must  analyze  I(D.  t)  and  I(t  -  Z   ).   I(D^  t) 
is  the  envelope  of  the  received  signal  (time  referenced  to  the  replica 
signal),  and  is  given  by 


X(D^  =  0/2  {K^t-iTo' ) - i[oAl -CUiTol}  ,  e-2 


Then  define 

T*    =  V/D  .  r_3 
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Therefore,  equation  E-2  becomes 


110*0  =  Otlx^iU-iTO-ltt-tU^T']]  .         M 


The  envelope  of  the  replica  signal  is  given  by 


M 


Kt-*)  =  o<  2oCi{iii-«r-iTp^  -  iii-T-a^U} .  e-4 


Figure  E-l  is  a  set  of  waveforms  of  the  received  and  replica  signals' 
envelopes »   We  have  chosen,  in  Figure  E-l,  a  code  of  length  6,  with  the 
code  values  of  4+l»  ~1»  +1»  +1»  -1»  +lj  as  an  example.   From  this  example 
we  will  obtain  an  equation  that  will  apply,  in  general,  to  codes  of  any 
length.   Note  that  the  amplitude  of  the  received  signal's  envelope  in 
Figure  E-l  is  less  than  the  amplitude  of  the  replica  signal's  envelope. 
This  is  done  deliberately  so  that  both  waveforms  may  be  seen  in  detail; 
however,  bear  in  mind  that  both  envelopes  are,  in  actuality,  of  the  same 
amplitude 

For  positive  doppler  we  recall  that  the  received  signal's  waveform 
will  be  compressed.   Figure  E-l  is  a  series  of  waveforms  of  I(D^  t)  and 
I(t),  in  which  we  hold  I(D^  t)  fixed  and  shift  I(t  -  T  )  to  the  right, 
corresponding  to  positive  range  deviation,  T  .   In  doing  so,  we  see  that 
the  limits  of  integration  (as  a  general  rule)  change  for  T  in  increments 
of  AT,  where 

At  =  1td1  -  r\     .  e-5 

Next,  we  would  like  to  express  the  interval  of  integration  (and  hence 
the  limits  of  integration)  in  the  form  of  a  series.   Let  us  now  look  for 
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symmetries  in  Figure  E-l.   Note  that  for  OC  X  <  4  AT,  the  limits  of 
integration  over  each  interval,  in  addition  to  the  code  values  of  the 
received  and  replica  signals,  remain  constant.   We  also  see  that  this 
occurs  again  for  9  ATC    X   «(13  AT0   However,  we  have  lost  the  duration 
of  one  code  symbol  of  the  replica  signal.   Similarly,  for  19^1^^^  23  AT, 
we  see  that  we  have  lost  the  duration  of  two  code  symbols.   From  this 
analysis,  we  can  obtain  a  general  formula  to  express  the  limits  of  inte- 
gration and  the  code  values  in  the  interval  of  integration.   hence  for  T 
in  the  interval 


T'   ,  r   ,  T'   ...  ,..   ,Trv 


l^Sr  ^-Sf  i-S-"^1-^^  >  --1.VI. 


AT  ~  AT 


AT 


where 


N  = 

Xi  = 

D.  - 

l 

D1  = 


D, 


TDi 


length  of  the  code 

Zl  depending  on  the  ith  value  of  the  code  symbol 

replica  signal's  doppler  factor 

received  signal's  doppler  factor 

differential  doppler  factor  =  D'/Di 

T/Dt 


TDi/D 


T' 
AT  -■   \TDi  ■  V\      , 
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The  wideband  single  code  ambiguity  function  becomes 
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in  which  the  definitions  in  equation  E-6  apply  and 

p  =  smallest  integer  >  \lt/"r'\ 

0(  =  normalization  constant  =  \  2±/\\T 


To  complete  the  analysis,  we  would  like  to  ohtain  the  code  values 
and  limits  of  integration  for  the  intervals,  4  AUT  <  8AT,  14  AT4/TO-8  AT, 
and  so  on.   As  a  general  rule,  the  interval  of  range  deviation  becomes 


T--](U4U-Blk  tX.    L    w\lL     '  m  -1,2,-",N. 

at  &t      at  ~*      &t 

We  see  that  we  may  express  the  ambiguity  function  for  this  interval  as 
the  sum  of  five  terms,  i„e8f 
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JJ-2> 


(JL+2.)T 


ic^4. 
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where  the  definitions  in  equations  E-6  and  E-8  apply  and 


n  = 


AT 


^ 


n=^i 
AT 


q  =  smallest  integer 

1  =  q  modulo  (n1   ) 

k  =  smallest  integer  >  \  7  A^..  \ 
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APPENDIX  F 
PROPERTIES  OF  COMPLEMENTARY  SEQUENCES  3,  15 

This  Appendix  is  a  compilation  of  some  of  the  properties  of  comple*- 
tnentary  sequences,   At  the  end  of  this  appendix  is  a  list  of  the  known 
kernels,  of  which  there  are  four, 

1.  Definition:   A  pair  of  hinary  sequences  of  equal  length,  with 
the  number  of  like  pairs  of  one  sequence,  A,  equal  to  the  number  of 
unlike  pairs  of  the  other  sequence,  B,  for  each  possible  spacing  except 
at  zero  spacing,  where  the  pairs  of  both  sequences  are  like,  is  said  to 
be  a  complementary  sequence,, 

2.  Given  sequences  A  and  B  each  of  length  n.   Then  complementing 
the  code,  i.e.,  complementing  each  element  in  the  code,  of  A  or  B  or 
both,  results  in  a  complementary  code. 

3.  Definition:   Time  inversing  of  a  code  is  to  interchange  the  1st 
and  last  bit,  the  second  and  next  to  the  last  bit,  and  so  on. 

4.  Time  inversing  the  A  code,  B  code  or  both  results  in  a  comple- 
mentary sequence. 

5.  Definition:   The  altering  of  a  code  is  a  transformation  where 
every  other  bit  of  both  codes  is  complemented. 

6.  The  result  of  altering  a  complementary  sequence  is  again  a 
complementary  sequence 0 

7.  Definition:   A  kernel  is  a  basic  length  code  which  cannot  be 
decomposed  into  shorter  length  codes  by  an  inversion  of  the  standard 
generating  methods. 

8.  Definition:   Complementary  sequences  which  are  not  kernels  are 
called  composite  complementary  sequences. 
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is  a  complementary  sequence  pair,  then 

form  a  complementary  sequence  pair. 

IUj      1  r   A  =   ''lijoctiocoo  »a 

!->      ni,.eoocoo,D 

is  a  complementary  sequence  pair,  then 

D  =  al»bl»a2»b2»  ° an»bn 

form  a  complementary  sequence  pair., 

11.   The  following  is  a  list  of  the  known  kernels: 

a.  length  2 

A  =  {-1,-lt 
B  =  {-1,+lJ 

b.  length  10 

Ax  =  {-l,+l,+l,-l,+l,-l,+l,+l,+l,-l} 
Bx  =  {-1,  +1,+1,+1,+1,+1,  +1,-1,-1,  +1} 

A2  =  £+l,-l,+l,-l,+l,+l,+l, +1,-1,-1} 
B2  =  {+1,+1,+1,+1,-1,+1,+1,-1,-1,+1} 

c.  length  26 

A  =  {+1,-1,+1,+1,-1,-1,+1, -1,-1, -1,-1, +1,-1, 
+1, -1,-1, -1,-1, +1, +1,-1, -1,-1, +1,-1, +1> 

B  =  {-1,+1,-1,-1,+1,+1,-1,+1,+1,+1,+1,-1,-1, 
-1, -1,-1, -1,-1, +1, +1,-1, -1,-1, +1,-1, +1> 


*T>±    indicates  the  complement  of  b^ 
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